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Abstract. Ultracold atomic quantum gases belong to the most exciting challenges 
of modern physics. Their theoretical description has drawn much from (semi-) 
classical field equations. These mean-field approximations are in general reliable 
for dilute gases in which the atoms collide only rarely with each other, and for 
situations where the gas is not too far from thermal equilibrium. With present-day 
technology it is, however, possible to drive and observe a system far away from 
equilibrium. Functional quantum field theory provides powerful tools to achieve 
both, analytical understanding and numerical computability, also in higher dimen- 
sions, of far-from-equilibrium quantum many-body dynamics. In the article, an 
outline of these approaches is given, including methods based on the two-particle 
irreducible effective action as well as on renormalisation-group theory. Their rela- 
tion to near-equilibrium kinetic theory is discussed, and the distinction between 
quantum and classical statistical fluctuations is shown to naturally emerge from 
the functional-integral description. Example applications to the evolution of an 
ultracold atomic Bose gas in one spatial dimension underline the power of the 
methods. 

1 Introduction 

Five millimeters of mere nothing separate the micrometer-scale ultracold cloud of sodium atoms 
from the glass walls of its surrounding vacuum cell in a typical experiment in the basement of 
maybe the reader's research institution. Five millimeters between the glass at 293 K and the 
Bose-Einstein condensed gas at a few nanokelvin. These eleven orders of magnitude in temper- 
ature are to be compared to the eight orders the temperature in a supernova is higher than 
that in our office. And the vacuum pressure of 10"^" Pa which is quantifying "mere nothing" 
is similar to the atmospheric pressure at the surface of the moon. 

BosG-Einstein condensation, a phenomenon predicted over 70 years ago by Albert Einstein 
[T] on the basis of Bose's new statistical formulation of a photon gas [5], has revolutionised 
atomic physics since its ground breaking experimental achievement in dilute alkali gases at 
JILA (Boulder) [3j, MIT [4], and Rice University [5]. By today, more than 80 groups world- 
wide dispose of techniques to produce ultracold and Bose-Einstein condensed gases, mostly 
of alkahs like ^^Rb, ^^Na, "^Li but also including metastable '^He, as well as "^^K, ^^Cr, ^^Rb, 
^^•^Cs, ^^''Yb. Another 100 groups need to be included when counting the wider range of ex- 
periments studying cold atomic gases, atom optics, trapping, cooling, and many more subjects 
[6] . The past decade has seen an exploding range of experiments studying many different prop- 
erties of such systems, thereby varying densities, atom numbers, dimensionality, interaction 
strengths, internal (electronic) state multiplicities, as well as character, geometry, shape, size 
and temporal behaviour of the external trapping potentials, both, between different runs and 
during the experiment's timeline. A range of review articles and monographs on theory and 
experiment can be consulted |7I8I9I10I11| . During the past five years, a growing number of 
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experiments with ultracold Fermi gases has been added to the spectrum of activities, see, e.g., 
Refs. |12|13|14|15|16j . For recent reviews see Ref. [17]. Today, ultracold fermions are regarded 
as a promising tool to design many-body quantum systems exhibiting many of the phenom- 
ena relevant in solid state systems, in particular in the context of superconductivity, and to 
explore these beyond the range of parameters realistic for degenerate electron gases. Moreover, 
new efforts aim at a cross-fertilisation between high-energy physics, in particular concerning 
the phase structure of quantumchromodynamics, and the physics of strongly correlated atomic 
gases. 

The quantum degeneracy of ultracold gases results from the (anti-)symmetrisation princi- 
ple for the wave function describing indistinguishable particles, together with the statistical 
behaviour of the many-body system. A non-interacting gas is characterised, for given mean 
energy and particle number, by the occupation numbers of the single-particle eigenmodes. 
Bose-Einstein condensation emerges as the macroscopic occupation of a single mode, at zero 
as well as finite temperatures |18) . Although interactions are not required for the existence of 
this quantum degeneracy, they are, in a realistic physical system, indispensable for reaching the 
respective degenerate equilibrated state. In fact, in experiment, the key last step to degeneracy 
is induced by evaporative cooling [7], where coUisional relaxation following a removal of the 
high-frequency tail of atoms restores an equilibrium distribution. 

Trapped atomic gases provide the unique possibility to tune both, the interaction strength 
between the particles and the external boundary conditions fixed by the trapping potential: 
External electromagnetic fields can be used to considerably vary, in particular near (Feshbach) 
zero-energy scattering resonances, the scattering length which quantifies the collisional interac- 
tions [19]. Thereby, laser light, combined with elaborate lensing technology, allows for almost 
arbitrary trapping geometries. Most of these tuning knobs can be turned so quickly as to excite 
many-body dynamics far away from a thermal or metastable equilibrium state. 

What do we mean by "far-from-equilibrium" dynamics? Consider a non-integrable many- 
body system. Far from equilibrium, in contrast to close to it, there is no longer the notion 
of a precisely defined spectrum of quasiparticle modes whose damping can be described on 
the grounds of a linear-response analysis. The latter generally relies on a perturbative expan- 
sion in some small parameter, and analytic relations between the fluctuation and response 
functions reflect the principle of detailed balance in Boltzmann's kinetic theory [20l21j . In far- 
from-equilibrium or, as it is also often called, nonequilibrium time evolution there is no such 
fluctuation-dissipation relation, while microcausality and microreversibility are still conserved. 
When studying the transition from far- from- to near-equilibrium dynamics one of the key ques- 
tions is how the known near-equilibrium features are recovered during the time evolution of 
the system, given specific interaction properties. An interesting observation gives rise to the 
distinction between short- and long-time evolution after a sudden quench of some boundary con- 
ditions which produces an initial state far from equilibrium. In particular long-time many-body 
evolution and equilibration are demanding and still largely unresolved problems. 

An important issue when studying nonequilibrium dynamics is the level of approximation 
on which interaction effects are taken into account. As pointed out above, interactions are 
inextricably linked with nonequilibrium dynamics. If they are weak, i.e., occur rarely which is 
the case in a dilute gas, low-order perturbative approximations in the diluteness parameter can 
provide a reliable description over a certain time. At large times, however, such perturbative 
descriptions are expected to break down. Moreover, if the interactions are sufficiently strong, 
the time at which this breakdown occurs can be shorter than the time at which near-equilibrium 
kinetic theory, which in general involves perturbative approximations, sets in to be valid. 

For systems with large, i.e., classical occupations of the kinematically relevant modes, quan- 
tum fluctuations typically play a minor role and, if classical fluctuations are relevant, Monte 
Carlo simulations are often the method of choice. However, if quantum fluctuations become im- 
portant, no such methods are at hand, and this is generically the case for long-time evolutions 
and dense, strongly interacting systems. (Note that recently, for a certain range of applications, 
stochastic simulation techniques have been studied |22l23l24l25j .) Functional quantum field 
theoretical techniques represent a powerful approach to such dynamics. Moreover, they provide 
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analytical insight and make numerical computations feasible, in particular in more than one 
spatial dimension. 

This article provides an introduction to functional quantum field theoretical methods to 
describe far-from-equilibrium many-body quantum dynamics. It is beyond its scope to rep- 
resent a full review and to provide a satisfactory account of the relevant literature. In Sec- 
tion [21 we will define the relevant observables and give a brief summary of mean-field the- 
ory for the short-time and near-equilibrium evolution of a Bose gas. In Section [3l we shall 
then discuss the nonequilibrium two-particle irreducible (2PI) effective action |26I27I28] in the 
nonperturbative approximation introduced in Refs. j29l30| . This nonequilibrium approach has 
been developed and extensively applied, in the context of relativistic quantum field theories, 
as reported in Refs. ^29 30 31 ■32 33 34 35136137] . Their extension to abelian and non-abelian 
gauge theories |38|39|40|41j is the subject of recent and ongoing research. We will, furthermore, 
summarise a new functional renormalisation-group approach to far-from-equilibrium dynamics 
introduced in [42]. For applications, the focus will be set on interacting ultracold atomic Bose 
gases |43|44|45|46| and an overview be given of the results first presented in Refs. [46|47|48|49j . 
Example applications to the long-time evolution of an interacting Bose gas will be discussed. 
The relation to near-equilibrium evolution is described in Section [4] Quantum statistical fluc- 
tuations and their distinction from their classical counterparts are the subject of Section [3 A 
summary will be given in Section [51 

2 Mean-field dynamics of Bose-Einstein condensates 

In this section we give an introduction to mean-field theory of time-evolving ultracold Bose 
gases. Mean-field theory is generically valid for the description of near-equilibrium dynamics 
and has been applied successfully to understand and predict an enormous variety of experi- 
mental observations, see, e.g. |10lll| . Discussing it allows us to lay the foundations for the later 
development of far-from-equilibrium dynamics and dynamics of strongly interacting systems. 
We will first define the observables of interest in the context of most experiments and then 
give a concise introduction to the Gross-Pitaevskii and Hartree-Foek-Bogoliubov theories of 
Bose-Einstein condensates near their ground-state configuration. We close the section with an 
outlook beyond mean-field theory, focusing on atomic gases near a Feshbach resonance as well 
as on an exact method to calculate the nonequilibrium dynamics of an interacting Bose gas in 
one spatial dimension. 

2.1 Observables 

The application of statistical many-body theory to quantum degenerate states, e.g. of non- 
interacting Bose-Einstein condensates (BEG), aims at the occupation number distribution over 
the available energy eigenstates. Quantum states, however, contain additional information which 
manifests itself in the phase of the wave function. This phase gives rise to coherence and inter- 
ference phenomena which, in a number of experiments with degenerate atomic gases, has been 
made visible in a macroscopic manner, e.g., in the experiments reported in Refs. [50l51l52l53l54j . 
Let us discuss in more detail the observables required to describe such properties of the many- 
body system. 

Large occupation numbers, together with the phase of the single-particle wave function, lead 
to an approximate description of coherent matter in terms of a scalar complex field 4>{x). Here 
and in the following time and space variables are included in the four- vector x = (t, x) = (xq, x). 
In view of the statistical properties of BEG, (j) forms an order parameter. With respect to 
coherence properties as well as large occupation numbers, the matter-wave field 4>{x) exists in 
full analogy with, e.g.. Maxwell's classical electromagnetic field Ffj_^{x). Recall that an ideal, 
single-mode laser can be described by a coherent state which characterises it as a coherent 
superposition of photon number states. 
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The expansion coefBcicnts are chosen such that \a) is an eigcnstate of the Fock annihilation 
operator a with eigenvahie a. For photons propagating in vacuum, the time evolution of the 
coherent state reads |Q;(t)) = |a(0) exp(— icji)), u being the frequency of the mode. This corre- 
sponds, in spatial representation, to a classical oscillation of the wave packet in the oscillator 
potential muJ^X^ /2. Hence, the expectation value of the electric field operator E oc £{a + a^) 
performs the harmonic oscillations of the macroscopic classical field. 

For an ideal-gas BEC, the field can be written as 4>{x) = ^/Nip{x), where ^{x) has the 
functional form of the quantum mechanical single-particle wave function whose macroscopic 
mean occupation number is A^. Hence, while |</'(a;)P describes the particle density at x, its 
phase gives rise to the same interference phenomena as the wave function ip does for single 
particles. The complex matter-wave field describes the spatial and temporal density and 
phase distributions of the coherent cloud of atoms in a particular internal state. This analogy 
with photons gives rise to the notion of an atom laser in experiments where, e.g., a coherent 
beam of particles is coupled out from a trapped BEC [51|52|53] . 

A remark is in order: Only compact systems are experimentally relevant. According to 
the above picture, is the expectation value of the complex non-relativistic quantum field 
operator obeying the bosonic equal-time commutation relations [^(x, t), <?'''(x', i)] = (5(x — x'). 
In finite, closed, non-relativistic gatherings of atoms, the total number of atoms is a conserved 
quantity, i.e., the expectation value of this field operator with respect to the reduced density 
matrix describing any subsystem necessarily vanishes, (j> = (^) = 0. Therefore, a coherent state 
can not describe such a system. This is closely related to the fact that an isolated system does 
by definition not interfere with any other system such that its total phase can not be measured. 
Phase can only be measured by means of interference effects, and hence only the relative phase 
between different (sub)systems is a meaningful quantity. 

This leads to the concept of the phase coherence of a BEC. This coherence manifests itself 
in the off-diagonal elements of the reduced single-particle density matrix, i.e., the single-time 
two-point correlation function 

n(x,y;i) = (.^*(x,t)<?(y,i)). (2) 

The local particle number density is given by the diagonal elements, n(x; t) ~ n(x, x; t). For an 
infinite uniform ideal-gas BEC, one finds that the first-order coherence function derived from 
the off-diagonal elements, n^^'>{s) = rt(R -I- s/2, R — s/2) const. ^ for s —> oo. Since the 
momentum distribution of particles is given by the spatial Fourier transform of n(x, y;<) with 
respect to the relative coordinate s = x — y, this asymptotic off-diagonal finiteness implies a 
macroscopic occupation of the zero mode, i.e., Bose-Einstein condensation. It is termed Off- 
Diagonal Long-Range Order (ODLRO) and allows the asymptotic factorisation of the single- 
particle density matrix, 

n(x,y;t) ^H^,my,t). (3) 

ODLRO was introduced by Penrose and Onsager as a general criterion for interacting BEC 
[55| . It is applicable also in non- uniform systems alternative to the field expectation value 4>{x) 
and remains meaningful for number-conserving states. The asymptotic factorisation in Eq. ([3]) 
defines (f) as an in general complex order parameter for BEC. In finite systems, once "long-range" 
is quantified, it can be taken as a measure for local order. It clearly expresses the fixed phase 
relation between distant points in the BEC. In the thermodynamic limit (j) can be identified 
with the field expectation value (j). 

In summary, the field expectation value (/) is a useful concept for most cases where a system 
involves macroscopic occupations TV 3> 1. When comparing a formulation in terms of (j) with 
that based on a density matrix with fixed total particle number differences in the observables 
are suppressed at least with a factor 1/N. 

In the following, the classical field 

<j>ix) = i'Pix)) ^ Tr[p{to)<Pix)] (4) 

will be regarded as the non- vanishing expectation value of the field operator evaluated at some 
time t with respect to the density operator p at some initial time to, where the time dependence 
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of <P is implied to include the time evolution from to to t. Besides this one-point function we 
will use, in the following, the time-ordered connected Green function 

Gab{x, y) = {T<Pa{x)'^b{y))c = {T$a{x)My)) - Mx)My) = {TSa{x)Sb{y)), (5) 

where T denotes time ordering of the subsequent operators and <P — (l> — (p is the operator of 
fluctuations around the classical field (j). For the first, the field indices are chosen to number 
= ^ and ^2 = as independent components, {a,b} G {1,2}. The time ordering allows us 
to write G as 

Gab{x,y) = Fab{x,y) - -sgn(.To - yo)pab{x,y), (6) 

where F and p involve the anticommutator and commutator of the fields, respectively, 

Fab{x,y) = ^{{<Pa{x),<Pb{y)})c, Pab{x,y) = i{[<Pa{x),'Pb{y)])c- (7) 

In the (P-<P^ -basis, one has Gn = {T<l>^(p)c, G21 = {T<P<P)c, etc., and G is related to the 
single-particle density matrix n = uq + n, riQ = and to the pair function m{x.,y;t) = 
(^(x, t)S{y, t)) as follows 

fn{x,y;t) = F2i{x,y)\^j:g=yg=f (8) 

As will be discussed in more detail in Sect. IH F is called the statistical correlation function, 
containing, near equilibrium, information about the occupation of the available states, while 
the spectral correlation function p provides the frequencies and widths of these states. Near 
equilibrium, F and p are linked by a fluctuation-dissipation relation, fh quantifies the amount 
of pair correlations in the system, as, e.g., the number of atoms bound in pairs or, for fermions, 
the number of long-range correlated (Cooper) pairs, see Sect. I2.4TT1 

Connected n-point functions with n > 3 contain information about higher-order correlations 
and are defined analogously. In this section we will concentrate on n < 2. 



2.2 The Gross-Pitaevskii equation 

This subsection intends to give a concise summary of the non-relativistic classical equation of 
motion for the field expectation value 0(a;) which represents the leading-order approximation 
to the full quantum dynamics in the case that 7^ 0. It was first studied by Gross [56j and 
Pitaevskii [57j in the context of vortices in superfluid helium. The equation for a field describ- 
ing an ideal-gas BEG according to our above discussion has the same form as the Schrodinger 
equation for the single-particle wave function ip and is a special case of the Gross-Pitaevskii 
equation (GPE) which, in addition, includes the leading-order effects of interactions between 
the particles. 

We define the system to be studied by the many-body Hamiltonian for a single-species of 
non-relativistic particles, 

H = y"dx<pt(x)i7iB(x)<P(x) + iy"dxdy<pt(x)<|)t(y)y(x_y)<^(y)<?(x), (9) 

with the Hamiltonian for a single particle exposed to an external, e.g., trapping potential VJ^xtj 

+ 2 

iJiB(x) = V^-f 14xt(x), (10) 
Zm 

and the two-body potential V{x — y) describing the interactions between the particles. The 
interactions are assumed to depend only on the relative coordinate between the colHsion partners 
as, e.g., in the Born-Oppenheimer approximation of atom-atom collisions. The GPE is then 
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obtained from the Liouville equation ihdt{^)t — {[^TH])t, where (•)( denotes the expectation 
value with respect to the density operator at time by approximating the 3-point function 
{(p^<P<P) as a product of field expectation values, {<P^(l»P) ^ ip*(j)(j>: 

ihdt4>{x) = [i/iB(x) + 4>{x). (11) 

Here, the contact-potential approximation has been chosen for V{r) = gS{r), with g ~ ATrh^a/m, 
where a is the s-wave scattering length. In Sect. 12. 4TT] below, we will discuss in more detail the 
scattering theory behind this parametrisation, and note here only, that this approximation ren- 
ders the GPE to be a low-energy effective theory: The typical length scale characterising the 
motion of atoms in a BEG, the thermal de Brogiie wave length At = /i/V^TrmfceT, is generally 
much larger than the scale determining the short-distance behaviour of the interatomic poten- 
tial. The atoms therefore only see an averaged interaction potential which at large internuclear 
distance manifests itself in a scattering phase shift. At low energies, given that the potential 
is sufficiently short ranged, only the s-wave scattering amplitude survives which tends to a 
constant related to the phase shift and equal to the minus the s-wave scattering length a at 
vanishing scattering momentum. The scattering length also quantifies the spatial extent of the 
highest excited dimer bound state close to the dissociation threshold and therefore the pair 
correlation length. Setting (7 = 0, Eq. (jlip has the form of the Schrodinger equation. 

We note that, in order to have Bose condensation in an ideal gas, \t needs to exceed the 
interatomic distance. In turn, for the GPE to be a good approximation, the interatomic distance 
must be much larger than the scattering length a. This means, the gas needs to be dilute, such 
that multiple scattering effects play no role. One says, the gas is weakly interacting. 

The quartic interaction term in the Hamiltonian shows, that the local energy density rises 
with increasing particle density if a is positive. On the other hand, if a is negative, the interaction 
part is lowered by increasing the local density of particles. Therefore, a positive scattering length 
is said to describe repulsive interactions while negative a corresponds to attractive ones. Note 
that despite this, also at positive scattering lengths the interatomic potential can support bound 
states and therefore exert attractive forces. Also the trap potential plays a crucial role as it, 
e.g., can stabilise a BEG of atoms with a < 0, since the kinetic part of the energy rises with 
increasing curvature of the field (f> at the density peak, see Eq. 

We briefly discuss stationary solutions of the GPE in the presence of a trapping potential, 
with the time dependence (j){t,x) ~ 0(0, x) exp(— z/it/7i), i.e., solutions of the time-independent 
GPE 

[i/iB(x) + .g|</)(x)P] 0(x) = ^</.(x). (12) 

For positive a, the dilute-gas BEG is often termed strongly interacting if the kinetic part can 
be neglected within the total energy. This is called the Thomas-Fermi (TF) regime, in which 
the density distribution obtained from (fT2|) reflects directly the shape of the potential: 

no(x) = |0(x)|2 = ^i^J^. (13) 
9 

Only at the edge of the cloud, where the density vanishes, the GP approximation breaks down. A 
BEG in a harmonic-oscillator trap with frequency Who is in the TF regime if the total number of 
particles times the scattering length is much larger than the oscillator length /ho — {'h/muji^oY/'^ ^ 
i.e., (3 = Na/lho ^ 1- The TF radius i?TF = P^^^hw characterising the extent of the cloud 
results from the requirement that N particles fit into the profile (fTS]) . Glearly, in the TF regime, 
the shape and size of the atom cloud is distinctly different from that of a nearly ideal gas, which 
is, to a good approximation, given by the modulus squared of the single-particle wave function. 

Gonsider, now, the GP dynamics disclosed by the non-linear field equation pT|) . Writing 
the field in terms of density and phase, '/'(a:) = \/ n{x) cxp{—iS{x)/h), one derives, by inserting 



^ "Strongly interacting" does here not imply that the gas requires a beyond-mean-field description, 
e.g., by being non-dilute in three or a geis with 7 > 1 in one spatial dimension, see Sect. 12.41 
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this into the GPE. the hydrodynamic equations 

5*71 + V • (nv) = 0, (14) 
m^t^r + V ( Vcxt + gn - ^^^^ V^^/^ + — j ^ (15) 

where the velocity field v is proportional to the gradient of the phase, v(a;) = {l/m)VS{x). 
Eq. (|14p is the continuity equation expressing local number conservation while Eq. (|15p rep- 
resents a quantum version of the Euler equation describing a frictionless fluid. The quantum 
contribution adding to the classical Euler equation reflects the zero-point fluctuations encoded 
in the term proportional to in Eq. (|15p : A curved mean-field profile is subject to a quantum 
pressure which aims at flattening the density distribution. We note that, since the velocity is 
a conservative or gradient field, the Euler equation describes, on a singly connected region of 
space, irrotational flow. As an example which exhibits in a nice way the consequences of this 
irrotationality consider the motion of a Bose-Einstein condensate which resembles that of a scis- 
sors mode first discussed in the context of nuclear physics. A condensed cloud trapped within 
an anisotropic, eUipsoidal potential which is excited by suddenly rotating the trap away from 
its prior position, starts oscillating around the new equilibrium orientation [58159160] . Although 
the oscillation of the density distribution resembles that of a rotation of the cloud, the velocity 
field shows that the flow pattern of particles is rather irrotational. 

Irrotationality of the flow is one signature of the superfluidity present in a system with Bose- 
Einstein condensation. A closely related and experimentally demonstrated property of such a 
system is the possibility of vortex formation. The Gross-Pitaevskii equation possesses non-linear 
solutions describing a circular flow around a singularity at which the density vanishes, i.e., 
the phase S (x) accumulates an integer multiple of 27r along one turn around the singular point 
(in 2 dimensions) or line (in 3D). Vortices can be excited in trapped BECs [61] using circular 
polarised laser beams [62] and have been observed to form Abrikosov lattice structures [63] as 
known from liquid Helium [64j . For a recent report on simulations see Ref. [65) . 

The GPE describes a colourful range of other nonlinear classical phenomena like solitons and 
nonlinear atom optics, phenomena which have been studied in many experiments and provide 
the frame of a research field in its own. See, e.g., Refs. |8lllj . With the advent of the formation 
of bosonic pairs in ultracold Fermi gases classical nonlinear dynamics can be studied in even 
more systems. 



2.3 Beyond the GPE: Hartree-Fock-Bogoliubov mean-field theory 

As discussed in the previous section, condensates exhibit superfluidity. The Gross-Pitaevskii 
equation for the order parameter field (t>{x) includes an Eulcr-like equation for a perfect fluid. 
Superfluidity, in turn, does not require a non-vanishing condensate order parameter as is known 
from the low-temperature physics of helium which is a non-dilute and therefore strongly in- 
teracting system. To describe BEG away from zero temperature and vanishing interactions, as 
well as away from thermal equilibrium, fluctuations need to be considered beyond the Gross- 
Pitaevskii approximation. In Section [2. II we identifled 4'{x) as a practically suitable measure for 
the asymptotic off-diagonal long-range order contained in the full two-point correlation function. 
Describing a Bose gas beyond the GP approximation requires additional information about the 
more local properties and dynamics of the two-point function, i.e., about the connected corre- 
lation function G{x, y) as well as about the back-reaction of this onto the evolution of (p. In the 
GPE the terms accounting for this back-reaction were neglected when approximating {<P^<P<P) 
as a product of field expectation values. 



2.3.1 Time-dependent HFB equations 



The dynamic equation for G can be derived, as before, from the Liouville equation ihdt{0)t ~ 
{[0,H])t, where O is now to be replaced by the respective products of field operators 'l>a{x)- 
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Fig. 1. Schematic representation of the different scattering processes contributing to the dynamics 
of a homogeneous gas within the Hartree-Fock-Bogohubov (HFB) mean-field approximation. The left 
panel illustrates the Gross-Pitaevskii (GP) approximation, which includes scattering between atoms in 
the condensate (p = 0) mode. The right panel distinguishes between the three possible channels for 
elastic scattering between a condensate and an excited atom, the direct (Hartree), exchange (Fock), 
and pair production (Bogoliubov) processes. In an operator language, these processes are described by 
the respective vertex operator contributions to the interaction Hamiltonian quoted above the diagrams. 
All vertex operators are at most quadratic in excited-mode operators flp (p 7^ 0) , such that the resulting 
Hamiltonian can be diagonalised, and it describes an effectively free system. 



Rewriting all correlation functions appearing in the equations in terms of their cumulants, i.e., 
in terms of connected Green functions and neglecting all cumulants of order three and higher 
one obtains the equations 

[ihdt - HiB (x)] 0(x, t) = g [(^(x, if + m(x, x; i)) 0* (x, t) + 2n(x, x; t) 0(x, t)] , (16) 
S^hdt ~ HiB (x) + HiB (y)] n(x, y;t) = g (^4){x, tf + 77i(x, x; t)^ m* (x, y; t) 

+ 2(0*(x,t)0(x,t)+n(x,x;t))n(x,y;t)] - .g[x^y]*, (17) 



[iWt - i?iB(x) - i?iB(y)] m(x,y;i) 



(x, tf + mix, x; t) ] n{y, x; t) 



H-2 (4>* (x, i)0(x, t) + 7i(x, x; i)) m(x, y; i) + <?[x 



(18) 



where i?iB(x) = — ?i^V^/2m -I- ycxt(x) denotes the one-body Hamiltonian. The generalised 
GPE and the equations ([T7]) and for the connected propagator contain, through the 
normal and anomalous density matrices n and m, respectively, only the statistical two-point 
function F, cf. Eqs. ([6|), ([8]). The above equations are commonly termed the (time-dependent) 
Hartree-Foek-Bogoliubov (HFB) equations [66167168] which describe the mean-field dynamics 
beyond the GPE in leading order in the coupling g. They form a closed system of partial 
differential equations which describe the coupled dynamics of the condensate and noncondensate 
components of an ultracold Bose gas. The set of equations preserves important conservation laws 
such as the total number of particles and energy. The exchange between the condensate and the 
noncondensed fractions is caused by the elastic direct and exchange collision processes between 
a condensate atom and an excited atom, as well as pair excitations out of the condensate, see 
Fig. [II 



2.3.2 Linearised HFB equations 



Let us finally consider the case that Eqs. ([T6|) ~([T8 |) describe small-amplitude oscillations around 
their stationary solutions. This approximation is, at first sight, irrelevant for the later discussion 
of far-from-equilibrium dynamics. We will discuss it here since it has widely been used for 
ultracold gases, and in order to point to the difference between the so called Landau and 
Beliaev damping processes and the collisional damping we will discuss in the later sections. 
See, e.g., Refs. [69|11| for a concise summary of the procedure outlined in the following. 
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On the right-hand side of Eq. (|16p . only diagonal elements of the normal and anomalous 
density matrices appear. We can therefore focus on the closed set of equations for <f){x) and the 
diagonal elements n(x, t) = n(x, x, t) and m(x, t) = to(x, x, t). To study small-amplitude devia- 
tions one linearises the generalised GPE (|16p as well as Eqs. (fT7|) and (fT8|) in small displacements 
from the equilibrium values of 0, n, and m, 

0(x,i) = (/)o(x) -h<5(/)(x,t), 
n(x, i) = n°(x) -I- <5n(x, t), 

m(x, = m°(x) (5m(x, f), (19) 

The time-independent part of the field expectation value is determined by the generalised 
stationary GPE 

fiM^) = [^ib(x) + g (7io(x) + 2H°(x) + m°(x))] M^), (20) 

with rio(x) = |(/)o(x)p. We work in the grand canonical ensemble, with the Hamiltonian replaced 
hy K = H — fiN which is equivalent to factor out a phase exp(— i/it/?;,) from (pQ in order to 
make it time-independent. 

In order to diagonalise the equations for the stationary densities H'^(x) and m'^(x) one 
transforms the fluctuation operators to a quasiparticle basis by means of the Bogoliubov trans- 
formation _ 

1^ <Pi^,t) (u,{^) «*(x) ^ (a,{t) 

Here, Uj and are quasiparticle operators which satisfy the Bose commutation relations 
[ai{t),Q^-(t)\ = 5ij, provided the mode functions u^, Vi are subject to the normalisation condi- 
tions / A'^x [u*(x)uj(x) - t;*(x)t;j(x)] = (5^ . 

Defining the normal and anomalous quasiparticle density matrices fij{t) — {al{t)aj{t)) and 
gij{t) — {ai{t)aj{t)), respectively, the Bogoliubov-deGennes eigenvalue problem 

/:(x) g [no(x) +toO(x)] \ / Uj(x)\ _ f Wj(x)\ , . 

-g [no(x) + mO(x)] -/:(x) J [v,{^) J ^ [v,{^) J ' ^^^> 

with C = i/iB(x) — /i + 2f;n(x), n(x) = ?t.o(x) + ^^"(x), fixes the quasiparticle amplitudes u,; and 
Vi, and yields a diagonal stationary part of the quasiparticle density matrix, f^j = fi^ij and a 
vanishing stationary anomalous quasiparticle density matrix g^^ = 0. 

The resulting linearised coupled equations for the time dependent variations (5(/)(x, t), Sfij (t) — 
fi]{t) - Sijff and 6g^j{t) g^j{t) read 

ihdtS(p{x,t) = (i?iB(x) - + 25ri(x)) 50(x,t) +.gno(x)(5(/)*(x,t) 

2,9(/)o(x)(5n(x,<) +g<?!)o(x)(5m(x,t), (23) 

tndtSf,,{t) = (e, - e,)<5/y(t) + 2g (/f - /]>) J d'^x 0o(x) [(<50(x, i) + 50* (x, t)) 

X (it,(x)u*(x) -f z;,(x)w*(x)) +(5</)(x,t)?;,(x)M*(x) +5(/)*(x,<)it,(x)u*(x)] (24) 

i/iat<55«,W - (e,+e^)%,(t) + 2.9(l + /f + /]>) J d^'x 0o(x) [((50(x, i) + (50*(x, t)) 

X (<(x)i;*(x) -f <(x)u*(x)) +<50(x,t)<(x)7.*(x) +<50*(x,i)<(x)«*(x)]. (25) 

The quantities /j* = (ajaj)o = [exp(/3ej) — 1]~^ are the equilibrium quasiparticle occupations, 
in terms of which the equilibrium non-condensate density is obtained as 

n^x) = ^[(|u,(x)p + |z;,(x)n/° + |«,(x)p]. (26) 



10 



T. Gasenzer 



The variations of the normal and anomalous particle densities expressed in terms of the quasi- 
particle density variations read 

,Sn(x,t) = ^ {K(x)^.,(x) + <(x)t;,(x)] 6f,,{t) + u,(x)t;,(x)%,(t) + <(x)^;*(x)<5g^(^)} , 

ij 

5m{^,t) - {2<(x)u,(x)J/,,(t) + u,{^)u,{^)5g,,{t) + <(x)«;(x)5g*.(t)} . (27) 

Note that Eqs. (PO]) . (|^^ - (P7)) form a closed set of equations of motion. 

Neglecting the cross terms coupling the condensate and noncondensate oscillations Scj) and 
5n, Sm, the equations for S(j) and S(f>* can be disentangled by a Bogoliubov-like transformation. 
The Bogoliubov frequencies are the resulting eigenvalues and therefore the frequencies of the 
BEC's elementary excitations. For a translationally invariant gas, they read 



7;— +9 no] -g rig 
2m 



1/2 

(28) 



where the Popov approximation = has been chosen j70|71j . This approximation ensures 
that the spectrum ([^ is gapless as required by the Hugenholtz-Pines [72] and Goldstone 
j73l74j theorems: At low momenta, the energy (j28[) is linear in the momentum p = |p|. In this 
limit the elementary oscillations are collective sound modes with dispersion ujp = CsP, where 
Cs = \J gno I m is the sound velocity. At high momenta, the Bogoliubov dispersion assumes the 
quadratic form of free particles, ujp = /2m + guQ. The approximations made in deriving the 
above equations are valid if the gas is weakly interacting, i.e., for a small diluteness parameter 
rj = na^ . 

One can show that, with a linear dispersion, energy and momentum conservation restrict the 
possibilities for the excitation of particle modes in a flowing BEG when encountering obstacles, 
e.g., atoms at a cavity wall (see, e.g., Ref. A linear dispersion implies a maximum critical 

velocity for frictionless flow, i.e., for superfluidity. A weakly interacting BEG therefore obeys 
Landau's criterion for superfluidity in the same way as superfluid ^He which, in addition, shows 
a pronounced roton minimum at finite wave vectors. Note that the Bogoliubov dispersion (j28[) 
is already obtained in the Bogoliubov approximation where the time evolution of the excited 
modes as well as the back action of the static excitation numbers on the condensate oscillation 
frequencies are neglected. 

We emphasise that the HFB equations ([M)) . ([^5]) are local in time and only involve single- 
time correlation functions. The full HFB equations (fT7|) . (fT8|l . however, also determine the 
off-diagonal time dependence of Gab{x, y). The fact that this does not feed back into the equa- 
tions P?)) . HH]) for the density matrices reflects that the HFB approximation does not account 
for direct scattering required for coUisional dissipation and thermalisation. Hence, the HFB 
approximation is expected to be valid in the coUisionless regime, where the mean free path is 
much larger than the scattering length. Note, however, that the linearised equations ((23l) - ((25)) . 
if the coupling of excitations of the condensate and noncondensed fractions account is taken into 
account, also describe one-to-two and two-to-one collision processes between the excitations, 
provided a BEG phase is present, i.e., (f>o /). These give rise to the so-called Landau and Beli- 
aev damping caused by the mixing of superfluid and normal fluid phases, see, e.g., |69|llj and 
Refs. cited therein. This damping is different in nature from the coUisional dissipation obtained 
beyond the HFB approximation of the dynamic equations and discussed further in Sect. [3?673l 



2.4 Beyond mean field 

The coUisionless regime discussed so far is easily left behind in present-day experiments. The 
preparation of ultracold atomic Bose and Fermi gases in various trapping environments al- 
lows to precisely study quantum many-body dynamics of strongly correlated systems, see, e.g., 
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Fig. 2. (Color online) Radial zero-energy scattering wave function rtl)o{r) = r limp^o V'p(i') (black solid 
line) at small internuclear distances r, for the square- well potential V{r) (drawn in blue). Shown are 
four different potential depths V'(O), resulting in different s-wave scattering lengths a. The scattering 
length is given by the intersection radius of the extrapolated wave function with the zero-energy axis, 
cf. Eq. (|29|l . (Red) horizontal lines in the potential wells indicate bound-state energy levels. 



Refs. |75l76l77l78j . In particular, techniques exploiting zero-energy (magnetic and photoasso- 
ciative, i.e., optical Feshbach) scattering resonances have helped to provide ultracold atomic 
gases with the importance they nowadays bear and the attraction they exert on physicists in 
most different areas of physics. Such techniques allow, by means of external electromagnetic 
fields, to tune the scattering length freely between large negative and large positive numbers. 
Special trapping configurations such as quasi one- and two-dimensional traps as well as optical 
lattices add to these possibilities and ask for descriptions beyond the mean-field level. 



2.4.1 Feshbach resonances 



Feshbach resonant scattering is most easily understood by realising that the value of the s- 
wave scattering length is directly related to the energy of the highest bound state the Born- 
Oppenheimer scattering potential supports below the zero-energy thresholdQ To illustrate this 
consider the simple case of a square-well potential V{r) which is non-zero only for r < ro 
and tends to infinity at r = 0, see Fig. [5] The radial zero-energy scattering wave function 
ripQ^r) = r limp^o ^'p(r) oscillates within the square well, with a frequency determined by the 
depth of the well, while its wave length outside the well is much larger, oc Hence, outside 



^ Overviews and references to most of the relevant literature concerning magnetic and optical Fesh- 
bach resonances and photoassociation can be found in Refs. [79180181 182I19I83I84I85] . 
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but close to the well, it is approximately linear in r. As can be seen in Fig. [21 the continuity 
conditions for the scattering wave function at the edge of the well imply that ipQ (r) crosses the 
r-axis at r = a. This can also be expressed in terms of the scattering amplitude approaching, 
for J) ^ 0, a constant, the s-wavc scattering length, linip^o fp{^) = ^o,' 



The square-well example shows clearly that, if the potential depth is changed, the energy Eh 
of the uppermost bound state shifts, and the scattering length goes through infinity when Eh 
crosses zero. We remark that the number of nodes within the well corresponds to the number of 
bound states left in the potential well. Very close to the resonance, the energy of the uppermost 
bound state with respect to threshold is proportional to the inverse of a^, 



The wave function of the uppermost bound state, as can be imagined from Fig. [21 is very similar 
to the wave function of the zero-energy scattering state for radii r considerably smaller than 
a. Close to r = a, however, since Ei, < 0, the bound-state wave function starts to differ from 
"00 (t") and approaches zero for larger r. As a consequence, the scattering length, if positive and 
larger than the extent of V{r), measures the spatial extent of the bound state, i.e., the size of 
the respective dimer molecules. Feshbach resonances have been exploited at large to produce 
degenerate molecular gases consisting of dimers of bosons, fermions, and of diatomic molecules 
of different species. See, e.g., Refs. |81|83I 84 85 86187] for reviews and further references on cold 
molecules. We remark that in these experiments, dimers could be identified at values of the 
scattering length on the order of one thousand Bohr radii [88] . These states are the largest and 
presumably most fragile molecules ever produced and measured in physics, see Fig. [3] 

In order to meet the conditions for a magnetic Feshbach resonance, the effective Born- 
Oppenheimer potential of the gas atoms is modified by means of external magnetic fields cou- 
pling to the magnetic moment. The atoms are usually trapped in a well-defined hyperfinc state, 
such that the field causes a Zeeman shift relative to the energy of atoms in different polarisation 
states. By applying an external magnetic field, different scattering channels, corresponding to 
different asymptotic hyperfinc states, can be shifted in energy relative to each other, see Fig. [U 
Intramolecular electromagnetic forces couple these potentials, with a strength depending on the 
internuclear distance as well as on the energies of bound states supported by the system. In this 
way, the effective scattering potential can be changed, and Feshbach resonances occur whenever 
a bound state of the coupled system crosses the energy of the asymptotically separated atom 
pair, see, e.g., Refs. |82|85] . We finally note, that photoassociation scattering is analogous to 
the magnetic Feshbach scattering described here. There, the coupling between the channels is 
provided by polarised laser light. In the most simple case of induced dipolc transitions, which 
require asymptotically P-wave closed channels, spontaneous decay of the closed-channel bound 
state results in a complex scattering length [90191192] . 

2.4.2 Ultracold gases near a Feshbach resonance 

In a many-body system, the description in terms of binary scattering becomes unreliable close 
to the Feshbach resonance. As follows from the above discussion, the scattering length can 
reach and exceed the mean atomic separation, and bound states are no longer binary but 
should be regarded as extended clusters involving a macroscopic number of particles. For short 
evolution times, theoretical comparisons with experimental results for molecule formation in 
condensates, see, e.g., Refs. (93,94 95 89 96]. indicate that the HFB dynamic equations discussed 
above can be applied even if the scattering length exceeds the mean interatomic spacing. In this 
experiment, performed in the group of Carl Wicman at JILA, Boulder, Bosc-Einstein condensed 
^^Rb atoms were, for the first time, observed to coherently bind to dimer molecules when the 




(29) 
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Fig. 3. Coupled-channel bound 
states of *"'Rb2 at magnetic field 
strengths of -Bevoivc ~ 16.0 mT 



and B„ 



= 15.55 mT. The dotted 



(dashed) curves indicate the closed 
({—2, —2} open) channel components. 
Note the extreme size at -Bmin- See 
text, Figs. \5\ M and Refs. |89I88| for 
more details. 



Fig. 4. (Color online) Coupled-channels picture of a mag- 
netic Feshbach resonance. The ultracold atoms collide, with 
an almost vanishing relative momentum, near the thresh- 
old of the background channel. Intramolecular forces couple 
them to bound states of the closed-channel potential Vd in 
which they are lacking the energy to separate to asymptotical 
distances. By tuning the Zeeman shift between the asymp- 
totic channel energies, a closed-channel bound state can be 
brought into resonance with the colliding atoms, causing a 
Feshbach-resonant increase of the s-wave scattering length a. 



scattering length a was tuned close to a Feshbach resonance [93] . The setup worked as a Ramsey 
interferometer: In the first step, the magnetic field B was tuned, for a few microseconds, close to 
Bq such that a{B) was on the order of 10'* Bohr radii qb- This constituted the first (coupling) 
Ramsey pulse. A longer evolution time ^evolve followed, during which a{B) was ramped back 
to a few hundred ab, and in the end, a second Ramsey pulse, identical to the first one, was 
applied, see Fig. O 

The length of the pulses were sufficiently short such that the HFB mean-field equations 
yield results quantitatively close to measurements [94I95I89I96] . Before the effects of multiple 
collisions and higher correlations become important a certain time after the quench a description 
within mean-field approximation remains valid. The calculations showed that the two Ramsey 
pulses coupled the colliding BEC atoms into the uppermost bound state of the two-channel 
system and coherently transferred atom pairs into molecules and vice versa. Due to the short 
pulse times, the coupling evolution resembled that of the fraction of a Rabi oscillation between 
two energy eigenstates. During the intermediate quasi free evolution, atoms and molecules could 
evolve such that a relative phase built up which was given by the binding energy relative to 
the free atoms multiplied by ^evolve- Depending on the value of this phase, the second Ramsey 
pulse lead to a further production of molecules or dissociated the previously formed dimers. 
As a result, sinusoidal oscillations of the remaining fraction of atoms at the end of the pulse 
sequence were observed as a function of iovoivc, see Fig. [6l The frequency of these oscillations 
precisely reproduced the expected binding energy of the Rubidium dimers [96] ■ 

However, discrepancies between theory and experiment remained in other cases, in particular 
for longer evolution times under strong interactions |97j , which may indicate that descriptions 
beyond mean field are required to interpret experimental data. We will discuss such methods 
in the following sections. 

Feshbach resonances have become, for experimenters, a versatile and convenient tool to 
control the coUisional interactions. The have, in particular, opened the door to the exploration 
of rich physics in parameter regimes never explored before. For example, ultracold atomic Fermi 
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Fig. 5. Scheme of a typical magnetic field 
pulse shape in the low density (no = 3.9 x 
lO'^^ cm"'') experiments in Ref. [93]. The min- 
imum magnetic field strength of the first and 
second pulse is Bmin ~ 15.55 mT. In the evo- 
lution period the field strength is chosen as 
Bevoivc = 16.0 mT. In the course of the ex- 
periments the evolution time fcvoivc as well 
as -Bovoivc were varied. The dashed line indi- 
cates the position of the resonance at Bo = 
15.49 mT. 



Fig. 6. The remaining fraction of condensate atoms, 
'ic(ifln), (solid line) together with the noncondensate 
fraction (dotted fine), and the total density of un- 
bound atoms (dashed line), as a function of the final 
time tfin, at the end of the magnetic-field pulse in 
Fig. [5l All densities are given relative to the initial 
density. The figure shows the results of simulations of 
the HFB dynamic equations in the form described in 
Ref. [89], for the experiment reported in Ref. [93]. See 
the main text for an outline of the experiment, and 
1891 for more details. 



gases can be manipulated such that they not only show the supcrconductor-like properties 
known from electron gases in solids but can cross from a BCS-like state containing Cooper 
pairs, over into a BEC of tightly bound molecules (For recent reviews cf. Refs. |98I99I17] ). As 
experiments are usually conducted in more than one spatial dimension and since the dynamics 
of Fermion gases can not be simulated by means of classical equations of motion, the functional 
field theoretical methods to be described in the following sections are expected to represent the 
most promising theoretical approach to the dynamics of strongly correlated Fermi gases beyond 
mean-field theory. 



2.4.3 Ultracold gases in lower dimensional traps and optical lattices 

One- an d two-dimen sional traps |11|76|77|100|101|102|103|104|105|106TT07] as well as optical 
lattices |108|109|110j allow to realise strongly correlated many-body states of atoms. In an 
optical lattice, strong effective interactions can be induced by suppressing the hopping between 
adjacent lattice sites and thus increasing the weight of the interaction relative to the kinetic 
energy [108llll| . This leads, in the limit of near-zero hopping or strong interactions, to a 
Mott-insulating state [75]. It is beyond the scope of this article to discuss in more detail the 
theory of ultracold gases trapped in such special configurations. Some remarks concerning one- 
dimensional (ID) gases, though, are in order, as we will focus on such a system when applying 
the functional field-theory methods in later sections. 

In special cases, the models describing ID gases }112lll3|114|115j allow to determine ex- 
act time-dependent solutions of the Schrodinger equation [1161117] providing insight beyond 
various approximations, which is particularly important in strongly correlated regimes. These 
ID systems are experimentally realized with atoms tightly confined in effectively ID waveg- 
uides [761771101110211031104110511061107] . where nonequilibrium dynamics is considerably af- 
fected by the kinematic restrictions of the geometry [77], while quantum effects are enhanced 
|118|119|120] . ID Bose gases are explored for various interaction strengths, from the Lieb-Liniger 
(LL) gas with finite coupHng [101I102I103I77] up to the so-called Tonks- Girardeau (TG) regime 
of "impenetrable-core" bosons |121|113|106|107|77] . The ID gas enters the TG regime if the 
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dimcnsionlcss interaction parameter 7 = giDmV/{h'^ N) is much larger than one. Here, giD is 
the coupling parameter of the one-dimensional gas, e.g., giu = 2h a/(nil\) for a cylindrical 
trap with transverse harmonic oscillator length /j^ In the Tonks-Girardcau limit 7 — > 00 
the atoms can no longer pass each other and behave in many respects like a one-dimensional 
ideal Fermi gas [TT] . 



2.4.4 Exact dynamics of an interacting ID Bose gas 



Most theoretical studies of the exact time-dependence address the Tonks-Girardeau (TG) 
regime |116|122|123|124|125|126|127|128| ). In this hmit, the complex many-body problem is con- 
siderably simplified due to the Fermi-Bose mapping property where dynamics follows a set of un- 
coupled single-particle Schrodinger equations [116j . A method for calculating the time-evolution 
of a LL gas with finite interaction strength has recently been discussed in Refs. [129ll30j . It has 
the potential to yield a valuable comparison to the results obtained from solving the dynamical 
field equations presented in the following sections. We therefore devote this subsection to a brief 
excursion and outline the method which generalises Lieb and Liniger's exact diagonalisation 
method [112j to time evolving A^-particle quantum mechanical wave functions. 

We consider the dynamics of N indistinguishable (5-interacting bosons in a ID geometry 
[112j . The Schrodinger equation for this system is usually written as 

idtipB = - g_^2 + 2cS{xi - Xj)'ipB, (31) 

i=l '* l<i<j<N 

where ipBi^i, ■ ■ ■ TXN,t) is the many-body wave function, and c quantifies the strength of the 
interaction, which is related to the dimensionless ID interaction parameter 7 = 2c introduced 
above. We do not impose any boundary conditions, i.e. , the x-space is infinite which corresponds 
to a number of interesting experimental situations where the gas is initially localized within a 
certain region of space and then allowed to freely evolve |123|124|125|126j . 

The idea is to construct exact solutions by differentiating a fully antisymmetric (fcrmionic) 
time-dependent wave function, which obeys the Schrodinger equation for a free Fermi gas |131j . 
The differential operator used for this depends on the interaction strength c and the number 
of particles. When c 00, the scheme reduces to Girardeau's time-dependent Fermi-Bose 
mapping |116| . valid for "impenetrable-core" bosons. 

Due to the Bose symmetry, it is sufficient to express the wave function ipB hi a single 
permutation sector of the configuration space, Ri : xi < X2 < ■ ■ ■ < xn ■ Within _Ri, ipB obeys 

N 

idt^B = -Y.9^^B/dxl (32) 
while interactions impose boundary conditions at the borders of i?i |112| : 

\ a a , ^B=0. (33) 

c \oxj-f-i ox 

This constraint creates a cusp in the many-body wave function when two particles touch, which 
should be present at any time during the dynamics. In the TG limit (i.e., when c 00) the cusp 
condition is ipBixi-, . . . ,Xj, Xj+i, . . . , xn, t)|a:j+i=2;j = |113I116| . which is trivially satisfied by 
an antisymmetric fcrmionic wave function ipp(xi, . . . , xn, t). Hence, ipB = within which 
is the famous Fermi-Bose mapping |113I116] . In many physically interesting cases, ipp can be 
constructed as a Slater determinant 

i^F{xi,...,XN,t) = (A^!)-^det[0„,(x,,O]^^,j=i. (34) 



^1/9 d 



16 



T. Gasenzer 



Since ipB — i^F within il)p must obey idi/jp/dt = ^ X^jLi ^^^V'F/^a;^, which imphcs that 
the (orthonormal) single-particle wave functions (xj , t) evolve according to 



Jdt = -d^cj)rn/dx^ 



(35) 



m = 1, . . . , N. Thus, in the TG limit, the complexity of the many-body dynamics is reduced 
to solving a simple set of uncoupled single-particle equations, while the interaction constraint 
([55)1 is satisfied by the Fermi-Bose construction. 

The simplicity and success of this idea motivates us to choose an ansatz which automatically 
satisfies constraint ([33|) for any finite c |131|115j . For this, define a differential operator 



O = 



n 



B 



(36) 



l<j<j"<Ar 



where Bij stands for 



Bi 



1 



1 



d 
dxi 



_d_ 

dxi 



It can be shown that the wave function 

V'B ^NcOipF (inside i?i), 



(37) 



(38) 



where A/'c is a normalization constant, obeys the cusp condition ([33| by construction |131|115| : 
Consider an auxiliary wave function 



'0Aux(a;i, ...,XN,t) = B-j+i^jOi/ji 



(39) 



where the primed operator Ojj_^_i ~ O/Bjj+i omits the factor Bjj+i as compared to O. The 
auxiliary function can be written as 



AUX 



1 



d 



dx 



_d_ 

dx. 



(40) 



It is straightforward to verify that the operator Bj^i jBjj^iO'j in front of "0^ is invariant un- 
der the exchange of Xj and Xj+i |115| . On the other hand, the fermionic wave function ijjp is anti- 
symmetric with respect to the interchange of xj and Xj+i. Thus, '0Aux(a;i, . . . , Xj,Xj+i, . . . , xjv, t) 
is antisymmetric with respect to the interchange of Xj and xj+i , which leads to 



V'AUx(2:i, • ■ • ,Xj,Xj+i, . . . ,XN,t)\2 



= 0. 



(41) 



This is fully equivalent to the cusp condition p3|) . Sj+i jV'-B Uj+ 1=3:3 = 0. Thus, the wave 
function ((38l) obeys constraint ((33)) by construction. 

In order to exactly describe the dynamics of LL gases, the wave function (|38|) should also 
obey Eq. ([32]) inside From the commutators [d'^/dxpO] = and [id/dt,0] = it follows 
that if ^/jp is given by Eq. and the (j),n{xj,t) obey Eq. ([55]) . then tpB obeys Eq. ([5^ . Note 

that for c 00, one recovers Girardeau's Fermi-Bose mapping |113|116j . i.e., O = 1. 

The above procedure can now be used to calculate the time evolution of a bosonic wave 
function which at time t = is given in terms of the operator Oc acting on some localised 
fermionic function. In order to calculate correlation functions, e.g., the single-particle density 
matrix ri, Eq. from this wave function, integrations over N ~2 variables must be performed. 
This can be achieved, at least for the diagonal elements of n, using Monte-Carlo techniques. 
We point out that the above method has the advantage that it does not rely on calculating 
the LL energy eigenfunctions of the system nor on the determination of the projection of the 
initial wave function on these eigenfunctions. The initial state is rather constructed from a 
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localised Fermi wave function. This approach sufBces for making benchmark comparisons with 
approximative methods as those studied in the following. 

We finally point out that for both, strongly correlated Bose and Fermi gases, methods in- 
spired by quantum information theory have recently been developed, with much effort and 
success, for the description of nonequilibrium dynamics. These time-dependent Density-Matrix 
Renormalisation-Group (tDMRG) methods draw from the empirical observation that the frac- 
tion of Hilbert space which becomes relevant during the time evolution of a generic many-body 
system scales only polynomially with the number of particle in contrast to the whole space 
which grows exponentially. For details on these techniques see e.g. Refs. [132113311341135) . 



3 Nonequilibrium quantum field theory 

The following section intends to give a concise introduction to the functional integral formula- 
tion of nonequilibrium many-body dynamics. The above brief discussion of mean-field dynamics 
has already given a taste of the formalism to be developed and made clear that the focus needs 
to be set on an efficient handling of the interactions which form the basis as well as the crucial 
difficulty of the whole theory. 

Assuming basic knowledge about Feynman path integrals we will sketch the functional 
integral formulation of real-time quantum field theory (QFT) and, in particular, introduce 
the one-particle irreducible (IPI) effective action. Eventually, to derive many-body dynamic 
equations beyond mean-field order which conserve crucial quantities like energy and particle 
number we will make use of an extended method which is based on the two-particle irreducible 
(2PI) effective action. We will consider initial-value problems in QFT which naturally lead to 
the concept of the Schwinger-Keldysh closed time path (CTP). The section will be closed with 
an application of the so derived equations of motion to describe equilibration of a Bose gas in 
one spatial dimension. 



3.1 Functional-integral approach 

We briefly recall the picture of quantum mechanics as provided by Feynman's path integral for- 
mulation. A common example is the time-evolution of a quantum mechanical state as exhibited 
by the transition amplitude from some initial state to a final state |tfin)- Given a classical 
action depending on the generalised coordinate ip(t) defined at each intermediate time step t 
between t-mi and tfin, the transition matrix element can be expressed in the path integral form 



defines the dynamical process in terms of the time interval to be considered and the Lagrangian 
L related to the Hamiltonian by a Legendre transform. The classical dynamics is determined 
through Hamilton's principle 




(42) 



where = 11*="* d,ip{t), and the initial 

valuefQ (^(tini) = {(filUni), V(^fin) = (<p|tfin)- 



and final states are reflected in the fixed boundary 
The classical action 




(43) 



6S[(t>] = 0, 



(44) 



^ We use the letter (p instead of x having in mind quantum mechanics as the 0+1 dimensional special 
case of field theory. 
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Fig. 7. (Color online) Classical vs. quantum 
mechanics. The classical path for given bound- 
ary conditions at fini and/or ta^ is shown as 
thick (red) line. The thin (black) paths would 
require, e.g., different initial values for ip, ip. 
In the microscopic world, the thin (black) 
paths add constructively to the path integral 
if their action S[ip] deviates less than h from 
the extremal value corresponding to the classi- 
cal path. Also tunneling processes as indicated 
by the thick (yellow) line would add construc- 
tively to the integral. 



which leads to the Euler-Lagrange, i.e., the sought dynamic equation for (p|f|For instance, given 
particular initial values for the position and velocity of the child on the slide shown in Fig. [7] 
at t = tini, this equation has the thick solid (red) path as solution. Different paths require in 
general different initial conditions to be imposed. 

On scales where quantum effects become relevant, the real world is somewhat more intricate. 
Fluctuations around the classical path as depicted by the thin (black) solid lines in Fig. [7] imply 
the action S[ip\ to deviate from its classical extremal value, and, only if this deviation is larger 
than h, the phase factor cxp{iS[ip]/h} suppresses the contributions of such paths to the integral 
through destructive interference. Qualitatively new effects are in order like the "quantum child" 
which can tunnel through the edge of the slide as along the (yellow) path in Fig. [71 

We generalise this path-integral formulation to QFT, where the coordinates (p become fields 
ip{x) defined over time and space. Moreover, we introduce external classical, i.e., non- fluctuating 
sources J{x) to turn the path integral into a generating functional for correlation functions, 
similarly as in the (grand) canonical partition function in equilibrium physics. This generating 
functional reads 

Z[J]^ j Vipe'^^^'f''^+i -^f^ (45) 

Here and in the following we shall use, if not explicitly stated otherwise, natural units, with 
h — \. We use the short-hand notation J Jkp — j^(i'^'^^xJ(x)tf{x) — J^^'" dxo J d'^xJ{x)(p{x), 
C = [iini,ifin]- For instancc, it allows the field expectation value (f) = (^) to be written as 

= Z-^ [ Vifi ip{x) e^'^^'^l (46) 

where W[J] = — i In Z[J] is the Schwinger functional. We introduce the quantum effective action 
r[(l)\ by demanding that the full quantum dynamics of the field expectation value (j) is given by 

In deriving the Euler-Lagrange equation the variation of the coordinate tp is usually taken to vanish 
at the boundaries of the time interval [tini,ten]. This procedure applies to systems with differential 
evolution equations of second order in time. For dynamic equations of first order in time, as the GPE, 
care needs to be taken when using the path integral for initial value problems, see, e.g., Ref. |136| 




SW[J] 



SJ{x) 



Ultracold gases far from equilibrium 



19 



Hamilton's principle applied to F, 

Sr[(l}] = 0. (47) 

This is equivalent to the functional measure containing a functional delta distribution which 
evaluates the fluctuating field to the expectation value cf) which is implicitly defined by Eq. (|46)) . 



Z[J]^ J V^S[ip-^]e''-^^'^^+^'^'^\ (48) 

with S[ip — (j>] = Yix^i'fii-'^) ~ '^(^))- Evaluating the functional integral shows that the IPI 
effective action defined as 

m - W[J] - I (49) 



satisfies Eq. (|48p . where it is implied that J, by inverting Eq. (j46|) . can be expressed in terms of 
4>. Eq. (|49p shows that the IPI effective action is the Legcndrc transformation of the Schwinger 
functional with respect to the classical source field J{x). While Eq. (|46l) defines the classical 
field as the derivative of W, the inverse expression (j47|) of J as a derivative of F constitutes 
the dynamic equation for (j) in the presence of the source J. 

Z is the generating functional for time-ordered n-point correlation functions, 



(Tc<?(xi)---^'(x„))= ^ ^"^["^^ 



Z[J] i5J{xi)---iSJ{xn) 



(50) 

j=o 



The Schwinger functional, in turn, generates the connected correlation functions or n-point 
cumulants. Hence, the connected propagator G defined in Eq. ([5]) is obtained as 



SJix)SJ{y) 



(51) 

j=o 



We briefly summarise the leading-order loop expansion of the IPI effective action r[(j)]. If 
the expansion of the classical action in powers of the field contains terms at most quadratic in 
and its derivatives, i.e., if there are no interactions present, the functional integral is Gaussian 
and therefore can be evaluated analytically. It is convenient to perform the shift (p ^ (p + (p 
and evaluate the Gaussian integrals over the fiuctuation fields (p. This yields the IPI effective 
action up to one-loop order. 



r("°°P)[0] = 5[(/.]-lTrlnGo, (52) 
where the classical two-point function Go is the inverse of 

*^o'^^'^) = ^50twl)- ^^^^ 

In the case that interactions, i.e., terms of cubic and higher order are present in the classical 
action S[(p], Eq. ([5^ represents the leading-order approximation to F if the corrections resulting 
from the interactions are sufficiently small. In this case the perturbative corrections can be 
calculated as a series of Feynman diagrams. We finally note that the second derivative of F[(j>] 
with respect to cf), 

' '^'^t'^^ -G-'{x,y)^G-\x,y)-S(x,y), (54) 



S(j){x)d(j){y) 



defines the proper self energy S to which only IPI diagrams contribute. 
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3.2 The two-particle irreducible (2PI) effective action 

The above steps set the framework of the functional approach to dynamics we choose in the 
foUowing. We will develop it slightly further and use the two-particle irreducible (2PI) effective 
action instead of the IPI one. The reason is that we are interested, as in the mean-field theory, 
in the time evolution of the two-point Greens function G{x,y). It is therefore desirable to obtain 
the respective dynamic equation in the same way as that for (f), as an Euler-Lagrange equation 
from Hamilton's principle applied to an action functional. The 2PI effective action which is a 
functional of both, (j> and G fulfills these requirements. Imagine it would be possible to calculate 
the IPI and 2PI effective actions exactly, e.g. by evaluating their loop expansions to infinite 
order. Results obtained from either of them would be identical. Their truncated expansions, 
however, are in general inequivalent. The loop expansion of the 2PI effective action can be 
understood as the sum of the diagrams contributing to the IPI effective action, reordered such 
that it can be expressed in terms of full propagators G which themselves can be represented as 
an infinite sum of diagrams involving only vertices and free propagators Gq. 

The 2PI effective action has been introduced to solid-state theory (there called (^-functional) 
in the sixties |26|27j and was later given its name in a relativistic quantum field theoretical for- 
mulation [28]. See also Refs. |137|138j . The nonperturbative 2PI 1/J\f approximation discussed 
in Sect. 13.6^ below was introduced in Refs. [29I30| . Detailed studied and applications to scalar 
relativistic as well as gauge theories can be found in Refs. [29|3 1 32 33 34 35 36 38 39 140(41] . to 
non-rclativistic systems, in particular, ultracold gases in Refs. [43„44.46.45i47.48.49j . 

To be specific, we consider a quantum field theory for a real A/'-component scalar field faix) 
(a = 1, with quartic interactions. Its classical action reads 



{x)ipaix)ipb{x)ipb{x), (55) 



with = J dxQ J d'^x. As we shall focus on a single-species ultracold Bose gas described by a 
complex scalar field, we choose Af ~ 2. The free inverse classical propagator then reads, in the 
basis where the field indices number the real and imaginary parts, </? = {tpi + iip2)/V2, 

iD~b{x,y) = jGnaf,(a;,y) =6{x-y) [-ial^dxa - i?iB (x)Jab] , (56) 

— 

see Eq. (|53|) . Here Hib{x) = —df /2m -I- V{x) denotes the single-particle Hamiltonian with 
interaction potential V{x) and the Pauli matrix in field-index space. Summation over double 
indices a, 6,c. . . — 1,2 is implied. With this, the action (|55]) corresponds to the Hamiltonian 
([5]), for a contact interaction potential V{r) = gS{r). The Bose field commutators in this basis 
read 

[■Pa{t,^),Mt,y)] = -^lbS{^- y)- (57) 

The definition of the 2PI effective action is based on a generating functional 

Z[J,R] = exp{iW[J,R]) = J ^y^e''^Sl^^+^^Mx)vA^)+^J^^v^(x)RMx.y)v,{y})_ (gg^ 

depending on the one- and two-point sources Ja{x) and Rab{x,y), respectively. The field ex- 
pectation value and the connected two-point function in the presence of sources are defined in 
accordance with the source-free case, 

?rrT-^"(-)' §^,-l(Mx)My) + G.bi..y)). (59) 

OJa[X) 5Rab[X,y) 2\ I 

The additional source term cx R can be seen as a quadratic (mass) term modifying S\^p\ — > 
S^\^p\ = S[(p] + J (pRip, such that the IPI effective action in the presence of R reads, to 1-loop 
order, cf. Eq. 

^i?,(iioop)j^j = + ^Trln [G^i(0) - iR] . (60) 
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Fig. 8. (Color online) (Left panel) Diagrammatic representation of the two- and three-loop diagrams 
contributing to the 2PI part 12 [4>, G\ of the 2PI effective action, cf. Eq. (|65|) . The bare vertices are drawn 
as black dots. Each such vertex is understood to represent a sum of the topologically different terms 
shown in the right panel. At each vertex, it is summed over double field indices and double space-time 
variables according to the respective diagram in the right panel. (Right panel) The representation of 
the bare vertex in terms of a black dot stands for a sum of the three topologically different connections 
of the four 'corners'. The black lines do not represent propagators and are only drawn in order to 
illustrate the different possible connections of propagators and/or external fields at the vertices. 



In order to arrive at the 2PI effective action one performs a further Legendre transform of 
with respect to the source i?, 



J XV oRab[x,y) 



r''[<P]- I cj,a{x)Rab{x,y)(j,k{y)-\TvRG. (61) 

J xy ^ 



Here we have used that 



^r'^m ^ 5W[J,R] ^ r / SW[,IR] _ \ SMz) ^ 6W[J,R] 
5Rab{x,y) 5Rab{x,y) SJc{z) " J 6Rab{x,y) 6Rab{x,y)' 



and its relation to G in Eq. ([59|l . From the double Legendre transform ([6T|) . p9|) . one directly 
finds the stationary conditions for (j) and G which, in the absence of sources, will provide us 
with the dynamic equations we are heading for: 

- -Ja(x)- / Rab(x,y)(l>b(y), j-„ , = --Rab(x,y). (63) 



S4>aix) Jy ' ' 5Gab{x,y) 

Before proceeding to these we derive the loop expansion of r[4>,G] which provides a way to 
obtain approximated expressions for the action and to perform practical calculations. 
Plugging Eq. ([501) ™to Eq. (|6ip yields the one-loop part 

r("°°P)[0,G] = S[c^] + ^Tr(lnG-i + GoiG) +const. (64) 
of the 2PI effective action 

r[<j,, G] = r("°°p) [0, G] + r2 [0, g], (65) 

where we have used Eq. ([5^ in leading-order, G^^ = Gg"^ —iR. Including all higher-loop terms 
in the "rest" /2 in Eq- ([65]) one finds, varying F with respect to G, that 

r,,(x,y;0,G) = 2*|^^|l^. (66) 
bGab\x,y) 



Since the self energy S is IPI, and since taking the derivative with respect to G corresponds to 
opening a propagator line, it follows that the rest term 7^2 must consist of 2PI diagrams only. 
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CTP C: 







+ branch 



- branch 



f 



Fig. 9. (Color online) Schwinger-Keldysh closed time path C. The green dots indicate the times xq 
and 3/0 for an example two-point function G{x,y), see text. The branches are drawn above and below 
the time axis only in order to make them separately visible. 



This forms the central result that the 2PI effective action is given, besides the terms (|64)) . by 
a series of all closed 2PI diagrams which can be formed from the full propagator G, the bare 
vertices defined by the classical action, and at most two external field insertions (/)■ 

The expansion of the 2PI part 12 [0, G] up to 3-loop order, for the classical action defined 
in Eq. ([55|l is shown in Fig. [51 We emphasise that, although the diagrams in this expansion are 
proportional to a power of the bare coupling g, truncations of the series can not be regarded as 
perturbative in g since the propagator G itself represents an expansion to infinite order in the 
coupling. The reason is that the stationarity condition for G, Eq. ([63]) yields a perturbatively 
truncated expression for the inverse of G. 

In order to arrive at a set of dynamic equations we need to discuss in more detail the 
implementation of the initial value problems we have in mind. 



3.3 Schwinger-Keldysh closed time path 

We assume that the many-body state is initially, i.e., at time t = to, given by some general 
(mixed) density matrix jo(to) The time evolution of the expectation value of an operator O is 
then given as 

{t\0\t) = Tr [p{to)U^{t, h)OU{t, to)] , (67) 

where U{t,t') ~ Texp{— i J* dt" H{t")/h} denotes the time evolution operator as obtained 
from the Hamiltonian H{t). 

The operators O relevant for us, i.e., the n-point correlation functions, are products involv- 
ing, in the Heisenberg picture, operators evaluated at different times. In the Schrodinger picture 
this implies additional time evolution operators between these factors. Consider, for instance, 
the two-time Green function, 

{Tc'Pa{x)My))c - Tr [p{to)Tc U\xo)^aix)U{xo) uHyo)My)U (yo)] - disc, (68) 

where U{t) = U{t^to) and the operators are time-ordered in a way which leaves the ordering 
within the products U^(l>U invariant. The disconnected part is denoted as 'disc.'. The product 
of different time evolution operators and field operators can be visualised by means of the closed 
time path as shown in Fig.[9l Starting at time to-, path sections leading to the maximum time 
appearing in the arguments of the field operators indicate time evolutions U. One generically 
chooses all times to lie on the -I- branch. However, different time orderings can be handled 
simultaneously by allowing times on the — branch as well and thereby doubling the range 
of possible times. Clearly, the two-point Green functions, with times evaluated on either or 
both of the two branches, are not completely independent from each other, and one aim of the 
discussion in later sections will be to clarify the dependencies. Here we only point out that the 
formalism to be developed naturally allows for two-time Green functions G{x, y) and therefore 
for Fourier transforms over their relative time xo — yo- These transforms, in turn, are functions of 
the frequency which, e.g., for a translationally invariant system, contain information about the 
spectral distribution of a particular momentum mode p. Beyond the mean-field approximation, 
collisions imply the redistribution of momentum between the particles. These scattering effects 
emerge naturally as finite widths in the spectral distribution around the dispersion peak at 
w(p). How these properties emerge from the dynamical theory to be developed is the topic of 
Section H 
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3.4 Gaussian initial states 



In the CTP formulation introduced in the previous section, the nonequihbrium generating 
functional for correlation functions can be split into a factor describing the initial conditions 
and one which contains all the ensuing quantum dynamics. For this we insert unit operators 
written in the coherent state basis: 

Z[J,i?,po] = / [d^+][d^o](^o+IPo|^o> V'^e^(^^'^^+f''^+^f^^^\ (69) 

The functional integrals over the fields at time io, [^''^^o'] = Ox '^¥'^(*0j x), evaluated on the + 
and — branches of the CTP include the initial density matrix. Their notion as limits of the 
dynamical functional integral imphes that the primed measure I?V na;o>to x d(f^ [x)dtf^ [x) 
excludes these initial-time fields. We emphasise that, due to causality, the CTP extends to 
the maximum time to be evaluated in a particular n-point function only. At later times, the 
sources can be set to zero such that the time evolution operators on the corresponding + 
and — branches cancel by unitarity. Note also that the CTP automatically arranges for the 
normalisation Z[0,0] = 1. 

The most general density matrix can be parametrised as 

{vt\P^W^)=^^'^'^''\ (70) 
with the normalisation factor M and fci'^] expanded in powers of the fields: 

oo „ ri 

^^i-^-<^'. mil 

Here, the coefficients a„ are non-zero only at the initial time to, at both ends of the CTP. It is 
now clear that Gaussian initial density matrices, for which = for n > 3, can be absorbed 
into the integrand of the dynamical integral in Eq. (j69p by a redefinition of the source fields 
J and R. The 2PI effective action approach then yields, once these modified sources are set to 
vanish, a closed set of dynamic equations for cj) and G. This set allows to specify initial values 
for these connected one- and two-point functions only, i.e., a Gaussian initial state implies that 
all higher-order connected n-point functions are assumed to vanish at i = toH 



3.5 Dynamic equations 



From the stationarity conditions (j63[) one obtains the dynamic equation for the field expectation 
value (f) which generalises the GPE to arbitrary order beyond the mean-field approximation, 
as well as a Schwingcr-Dyson-like equation which, upon multiplication by G has the form of a 
time evolution equation for the two-point Green function G, 

Go^lci^,z)Gcb{z,y) = SabSc{x - y) + / [Sac{x, z) + iRac{x, z)]Gcb{z,y), (72) 

J z 

with Sc{x — y) ~ 6c{xq — yo)6{x — y), which is a first-order differential equation in xq for 
the non-relativistic case with free propagator Using the definition ^ of G in terms of 

the statistical and spectral functions as well as the commutation relations ([57|) one derives the 



^ In the case that, at t = to, the nth-order connected correlation function {Tc$ai{to,^i) ■ ■ ■ 
<fa„ (io, x„))c is non-zero, with all mth-order functions, m > n, vanishing, a straightforward gener- 
alisation of the approach involving the nPI effective action is at hand [38] . 
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following set of dynamic equations: 

( - ^(^Ib^xo - 9Fab{x,x)^4)b{x) - {Hib{x) + ^[<t>c{x)(t)c{x) + Ft,c{x,x)\^4>a{x) 

\yEP^^{x,v-c^ = Q)c^b{y). (73) 



'to 

r'Xa ryo 

[^'^lc9xo+Macix)]Fcbix,y)^- dzX:P^{x,z;(l))Fcb{z,y)+ / dz Sj^^ix, z;(t))pcb{z,y), 

J to J to 



xo 



[i'^lAo + Mac{x)] pabix, y) ^ - dz S^^(x, z; (j))pcb{z, y), (74) 



yo 

t' t' 

Here we employ the notation dz = dzo / d'^z. The "mass" matrix M is defined as 

Mab{x) = 6ab Hib{x) + ^ (^M^) M^) + F^dx , x)'^ + .9 (a;)(/)fc (x) + i^ab (x, x)) (75) 

Moreover, the self energy has been decomposed into a part local in time and a nonlocal part 
written in terms of statistical and spectral components. 

Kb{x,y) = i:^"'(x)(5c(x -y) + ^ffc(x,y) - ^sgnc(xo - yo)S^bi^,y). (76) 

The local part, Z'^^''(x), has been absorbed into the mass matrix M while the non-local parts 
form the kernels for the memory integrals on the right-hand sides of the intcgro-diffcrcntial 
dynamic equations ([73|l . ([74|) . Note that only the two-loop, double-bubble term in / 2 contributes 

to 

With and UP set to zero, which is equivalent to truncating the loop expansion of the 
2PI effective action after the double-bubble diagram, one obtains a set of differential equations 
which are local in time and can be shown to be equivalent to the HFB dynamic equations 
introduced in Sect. 12.31 Note that the time derivatives in the equations for F and p only act 
on the respective first time variables xq. The corresponding equations in the second variable yo 
are obtained by symmetry considerations for F and p. The equations for different yo decouple, 
and only the diagonal functions with xq = yo appear in M. As initial conditions, (j) and F need 
to be specified at xo = yo = to- Note that p(to;x, y) is fixed by the commutation relations, 
Eq. dllD. 

Hence, in order to obtain Eqs. p7|) and (fT8|) one needs to combine the equations ([74]) . for 
Xo = yo with their counterparts with time derivatives acting on yo- The equation for p which 
is decoupled from all other equation is thereby shown to be consistent with the conservation 
of the Bose commutators. In summary, the HFB approximation is represented by the leading 
order (double-bubble) diagram contributing to /2. 

Setting, furthermore, F(to,x;to,y) = 0, -F remains identically zero at all times, and the 
equation for </> reduces to the Gross-Pitaevskii equation, cf. Eq. (fTT|) . 

The non-Markovian integral terms on the right hand sides of Eqs. ([73]) and (|74p open in 
a clear way the path to noncquilibrium dynamics beyond the mean-field approximation. The 
higher-order loop diagrams, with one line opened yield contributions to the self energy which de- 
scribe scattering processes involving the redistribution of particles between the different modes 
of the system. As the corresponding contributions to E^'P are non-local in time they natu- 
rally imply the evolution of a non-trivial relative time dependence of the two-point correlation 
functions and therefore finite-width spectral distributions as discussed above. 

Higher-order correlations formed in this way are not included explicitly in terms of higher- 
order n-point functions. One may, however, imagine these functions to obtain non-zero values 
implicitly, according to their own dynamic equations which then have been integrated formally 
and reinserted into the equations for the one- and two-point functions. 

Before we discuss in more detail the analytical implications of the non-Markovian terms 
and relate them to results obtained in the framework of kinetic theory we briefly discuss, in the 
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following sections, specific truncations of the 2PI effective action and, in a particular truncation, 
an application to a weakly interacting ultracold Bose gas evolving from a nonequilibrium initial 
state in one spatial dimension. 



3.6 Truncations of the 2PI effective action 

In order to practically solve the dynamic equations ([73|l and (|74|). details about the self energy 
S are required, and these are, in general, only available to a certain approximation. As discussed 
in Sect. [3?2l the natural expansion of /2 is in terms of 2PI closed loop diagrams involving only 
bare vertices and full propagators G. This expansion can be truncated at any order, e.g., of 
powers of g or the number of loops, without violating the most important conservation laws for 
energy and particle number which we briefly discuss in the following. 



3.6.1 Number conservation 



Particle number conservation is a consequence of the Noether theorem in conjunction with the 
invariancc of the theory under orthogonal transformations and can be seen as followsH The 
stationarity conditions (j63p can be combined to the equation 



'ab 



54>b[x) 



(5r[0, G] 

SGcb{y,x) 



Gcaiy,x) 



0, 



(77) 



with the elements CT2.a6 of the Pauli 2-matrix. From the specific expression (pS]) for the 2PI 
effective action follows that Eq. ([77]) is equivalent to the relation 



dxa-nix) - Vj(x) 



(f>a{x) 



^-Hntj^, G] 
S(j)bix) 



SFir^tjcl), G] 

SGcb{y,x) 



Gca{y,x) 



Here, 



n{x) = (j)aix)(j)a{x) +Gaa{x,x), 

j(x) = - [Mx)yMx) - Mx)yMx) + {rciS2ix)vSi{x) ^ Si{x)vS2ix))), 



m 



(78) 



(79) 
(80) 



are the total number and current densities, respectively. Clearly, particle number is conserved 
locally if n and j obey a continuity equation, i.e., if the right hand side of Eq. (|78p vanishes 
identically. We consider the specific structure of these terms: The interaction part of the 2PI 
effective action occurring therein is defined as 



,G]^n 



Mx)^DJ{x,y)My) - ^Tr [D-'G] 



xy 



The 2PI effective action is, like the underlying classical action S[ip] (|55p . a singlet under 0(2) 
rotations. It is parametrised by the fields <f>a and Gab, where the number of ^-fields has to be 
even in order to construct an 0(2)-singlet. From the fields (f>a alone one can construct only one 
independent invariant under 0(2) rotations, which can be taken as tr{(p(j>) = (jp' = 4'a't'a- AH 
functions of (p and G, which are singlets under 0(2), can be built from the irreducible, i.e., in 
field-index space not factorisablc, invariants [29|30j 



tr(G"), 



and 



tr(0</.G"), 



(81) 



Conservation of total particle number in a non-relativistic system corresponds to conservation of 
the difference of particles and antiparticles in a fully relativistic approach, i.e., to the conservation of 
charge. Neglecting the antiparticle sector of the Hilbert space in a non-relativistic system is equivalent 
to the constraint of a vanishing antiparticle number and, hence, amounts to the fact that the 0(2) 
symmetry of the Lagrangian ensures total number conservation. 
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with n = 1,2,.... As before, the trace tr(-) only applies to the field-component indices while 
there is no integration over space-time, e.g., tr(G^) = Gabix,y)Giic{y, z)Gcaiz,x). 

For contributions to /Int which contain only (jp or tr(G'"), the terms in square brackets in 
Eq. (|78p either vanish separately or are symmetric under the exchange of a and b. Moreover, 
if a term contains an invariant of the form tr((/)c/)G"), as, e.g., the contributions remaining in 
/],it from Tr{ G^^G}, the combination of the terms in square brackets in Eq. (j78p is symmetric 
under transposition in field index space. Hence, the total number density is conserved locally 
as a consequence of the 0(2) symmetry of the theory, and, more importantly, this is true for 
any set of approximate dynamic equations derived from a truncated but still 0(2)-symmetric 
effective action. Note, finally, that only terms in the action which contain mixed invariants 
tr((/)0G") induce exchange of particles between the condensate and the non-condensed fraction 
of the gas. 

3.6.2 Energy conservation 

Energy conservation follows from time translation invariance of cf., e.g., Ref. [139j . Consider 
the general translations in continuous space and time which vanish at the boundary, — *■ 
x'^+e'^{x), where e^(a;) is a time- and space-dependent infinitesimal 4- vector. The mean field and 
2-point functions transform, under these translations, to leading order in e, as 4>a{x) — > (t>aix) + 
£''{x)dl4)a{x), and Gab{x,y) Gab{x,y) + e''{x)d^Gab{x,y) + {y)dlGab{x,y).. respectively. 
Here, df, = d/dx'^, etc. One can show that under these transformations the variation of the 2PI 
effective action F can be written as F[(f>, G] F[4>, G] + SF[(t), G], with 

6F[^,G] = / T'^-'ix)d;e,{x). (82) 

Since, by virtue of the stationarity conditions (|63p . the variation SF vanishes for all solutions of 
the equations of motion for (jja and Gab, an integration by parts shows that T^'' is the conserved 
Noether current for the time-space-translations: 

SF[^, G]^- f e,{x) dlT'^'^ix) = 0. (83) 

J X 

F^^{x) is identified as the energy- momentum tensor, and the conservation law for total energy 
is expressed as d^T^^{x) = or dt J d'^xT^'^{t,x) = 0. Explicit expressions for the energy- 
momentum tensor have been calculated in Refs. |139I47| 

3.6.3 NLO 2PI 1/JV expansion 

Also the expansion of I2 in terms of 2PI loop diagrams can be resummed to obtain alternative 
non-perturbative approximation schemes. The most outstanding such scheme is the expansion 
in powers of the inverse number of field components Af applied to the dynamical 2PI effective 
action first by Berges and collaborators [29,30 39; and extensively studied since. To next-to- 
leading order (NLO) this resummation scheme can be understood as the replacement of certain 
vertices in a loop expansion by a bubble-resummed vertex |30I46| . The result of the scheme 
has also been recovered using a functional renormalisation group inspired approach |42| where 
it results as a truncation in orders of proper n-point functions combined with an s-channcl 
approximation of the equation for the proper four- vertex. 

In the following sections we will employ the NLO 1/Af expansion scheme to calculate the 
dynamics of a single-species ultracold Bosc gas described by a complex scalar field for which 
A/" = 2. In the context of a non-relativistic Bose gas, this approximation has been discussed in 
detail in Refs. [46|47l48|49j . In this scheme, the contribution F2[(f>, G] to the 2PI effective action 
involves a leading (LO) and next-to-leading order (NLO) part which can be diagrammatically 
represented as shown in Fig. [TOl While the leading-order contribution involves one diagram. 
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Fig. 10. (Color online) Diagrammatic representation of the leading order (LO) and next-to-leading 
order (NLO) contributions in the l/A/'-expansion, to the 2PI part /2[(/), G] of the 2PI effective action. 
The thick blue lines represent 2-point functions Gab{x,y), the red crosses field insertions (f>a{x), and 
the wiggly lines a single vertex channel of those shown in the right panel of Fig. |8] At each vertex, it 
is summed over double field indices a and integrated/summed over double time and space variables x. 

in NLO a chain of bubble diagrams is resummed. All of these diagrams are proportional to 
the same power of l/A/" since each vertex scales with l/A/", which is cancelled by the (blue) 
propagator loops which scale with J\f since they involve a summation over the field indices from 
1 to A/". Note that the Hartree-Fock-Bogoliubov (HFB) approximation is given by an action /2 
which involves and the first diagram of F^^'^ in Fig. [10] (b), see, e.g., Ref. |46j . 

From r2[(t>, G] = T^^^i^, G]+rf^'^[(j), G] we obtain, using Eq. the self energies Eab{x, y) = 
^abi^^y) - (i/2)sgnc(a;o - yo)E^Ax,y), with 
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AF{x,y) \Ap{x,y) 
-\Ap{x,y) AF{x,y) 



Fab{x,y) 
-^Pab{x,y) 



(84) 



where Ap^p(x,y) = lF,p{x,y) + Pp^p{x,y] If,p)- The resummation to NLO in \/M is expressed 
by the coupled integral equations for If,p [39] : 



lF{x,y)\ ^ g_ 
Ipix,y) J M 



F{x,yY -\p{x,yY 
2Fab{x,y)pab{x,y) 



If{x, z) 
Ip{x,z) 



dzl, 



ta 



'^^'^^ 1^ 2Fab{z,y)pab{z,y) ) 



2Fab{z,y)pab{z,y) 



(85) 



Here, = FabFab, etc. The functions Pf,p, which contribute to Af,p in the self energies 
and vanish ii (f>i = 0. read HSj 



2g ( n° 

PF{x,y;lF,p) = -jj:\^HF{x,y) + J dz [Hf{x, z)Ip{z,y) + If{x, z)Hp{z,y)] 

dz[Hp{x,z)lF{z,y) + Ip{x,z)HF{z,y)] - / dv dw Ip{x,v)HF{v,w)Ip{w,y) 

to ^ to Jto 

dv / dw Ip[x,v)Hp[v,w)lF[w,y) + I dv I 

to Jto Jto J Vo 



dwlFix,v)Hp{v,w)Ip{w,y)> [86) 



P. 



ix,y;lF,p) = -^{Hp{x,y) - I dz[Hp{x,z)Ip{z,y) + Ip{x,z)Hp{z,y)] 
■'^ Jyo 



xo pvo 

dv / dw I, 

yo Jyo 



p{x,v)Hp{v,w)Ip{w,y)Y 



(87) 
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wherein the functions Hp ^ are defined as 

HF{x,y) = -(j>a{x)Fab{x,y)(f>b{y), 

Hp{x,y) = -<f>a{x)pab{x,y)(f>b{y)- (88) 

The technical procedure to solve the above dynamic equations in every time step requires the 
determination of the functions I(x, y), before the actual propagation of the respective correlation 
functions. 



3.7 Functional renormalisation-group approach 

We close our introduction to noncquilibrium quantum field theory with a short description of 
an alternative approach which is based on functional renormalisation group (RG) techniques. 
For more details on this approach see Refs. [42|140j . Dynamic equations will be derived which 
are similar in structure to the equations obtained from the 2PI effective action. In a particular, 
"s-channel" approximation these equations are equivalent to the 2PI equations in NLO of the 
1/J\f approximation summarised above. The functional renormalisation group (RG) techniques 
we will employ have been introduced and used extensively in the framework of equilibrium 
quantum field theory of strongly correlated systems. See Refs. |141|142|143|144j . as well as 
|145j for non-equilibrium applications. 

For a given initial-state density matrix po{to), the renormalised finite quantum generating 
functional for time-dependent n-point correlation functions, 

Z[J-pd] =Tr[pDito)Tccxp{i J^^Ja{x)^aix)]], (89) 

(summation over double indices is implied) carries all the information of the quantum many- 
body evolution at times greater than the initial time tg. 7c denotes, as before, time-ordering 
along the Schwinger-Kcldysh closed time path (CTP). sec Sect. 13.31 and ^ = J^^dxo J d'^x. 
All connected Greens functions will be time-ordered along C. 



3.7.1 Functional flow equation 

The key idea of the approach to be described in the following is to first consider the generating 
functional for Green functions where all times are smaller than a maximum time t. This implies 
a time path C(r) which is closed at t = r, i.e., the maximum time in Fig. [9] is set to r. As a 
consequence, the generating functional Zr = ^c(t) has the source term 

rc(,)exp|z / Ja{x)^a{x)}. (90) 

^ Jx,C(r) J 

At r = toi this results in a trivial Zt^ where all information is stored in the initial density 
matrix poito)- From this initial condition Zr can be computed by means of the time evolution 
drZ-r for all times t > to. 

In the following, this evolution shall be derived by use of functional RG ideas. We note that 
Zt can be defined in terms of the full generating functional Z^o in Eq. ((89|) by suppressing the 
propagation for times greater than t. This suppression is achieved by 

Zr = cxp{-'-f ^^Rr^,,{x,y)-^\Z, (91) 

^JxyfiOJaix) 5Jb[y)i 

where the function is chosen such that it suppresses the fields, i.e., 5/dJa, for all times t > t. 
This does not fix i?^ in a unique way, and a simple choice is 

I oo for xn = yn > T, X = y, a — b 
otherwise 
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Fig. 11. (Color online) The cutoff 
function R^^ab{x, y) in the time plane 
{xo,yo] = {t,t'], t,t' > to- The func- 
tion vanishes everywhere except for 
t = t' > T where it tends to infinity 
and therefore implies a suppression of 
all fluctuations in the generating func- 
tional at times greater than r. 



Fig. 12. (color online) Diagrammatic representation of the 



general flow equations for J'r^'' [0], 



-.(2)r 



and r 



(4) 



0], for a (fi'^-theoiy. Open circles with a r denote drRr.ab, 
solid lines with (blue) filled circles are r- and, in general, (j)- 
dependent two-point functions GT,ab = i[rS^'' + Rt]~j^ ■ All 
other filled circles denote proper field-dependent n-vertices 



r. 



(n) 



3, 4, 6. P implies a sum corresponding to all 



permutations of its arguments. 



see Fig. [TT] 

We emphasise that the cutoff Rr in Eq. ((9T|) suppresses any time evolution at times greater 
than r. Correlation functions derived from vanish as soon as at least one of their time 
arguments is larger than t. Hence, the regularised generating functional Eq. (|9ip is equivalent 
to a generating functional with a closed time path C(t) leading from to to r and back to to- 
Note that the CTP automatically arranges for the normalisation of Zr- 

The restriction of the CTP to times to < t < t implies that the differential equation for Zr 
describing the flow of the generating functional, and therefore that of the correlation functions, 
encodes the full time evolution of the system. Analogously, the time evolution of connected 
correlation functions is derived from that of the Schwingcr functional Wr = ~i In Zt-. It is more 
convenient, however, to work with the effective action 



-^t[0;-Rt] = Wr[J;PD] - J Ja4>a ^ 2 j ' 



(93) 



Here, space-time arguments are suppressed, and 4>a{x) = 5Wr / 5Ja(x)\j=o is the usual field 
expectation value. From Eqs. ([OTI) and ([M)) we derive the Functional RG or flow equation for 
the r-dependent effective action. 



drRr^ab 



(94) 



- ab 



where Fr^^ = S"^rr/{6(j))". Again, space-time arguments are suppressed which appear in analogy 
to the field indices a, 6, see e.g. Ref. |144| . Eq. (|M|) is analogous to functional flow equations 
used extensively with regulators in momentum and/or frequency space to describe strongly 
correlated systems near equilibrium [141114211431144] . Its homogenous part relates to standard 
T-dependent renormahsation [144j . and has been studied e.g. in |146ll47j . 
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3.7.2 Flow equations for correlation functions 

To obtain a practically solvable set of dynamic equations, we derive the flow equation for the 

proper n-point Green function I^t"'' by taking the nth field derivative of Eq. ((94|) . Fig. [T2l shows 
a diagrammatic representation of the resulting equations for the r-dependent proper two- and 
four-point functions. To be more specific, we consider, in the following, the special case of an 
A/'-component scalar 0^ theory defined by the classical action ([55|) . 

Our goal is to derive the full time-evolution of -T*-"-* = in particular, of the connected 

two-point fimction G = = i[r^'^]-^. As a consequence of the above mentioned effective 

cut off of the CTP at times greater than r, it will be sufficient, for the time evolution up to 
< = r, to determine the functions I^r"^ and thus the propagator 

Gr,ab = i[r^^^ +Rr]ab- (95) 

For conciseness we here only discuss the case (f>a = 0, such that the action ((55|) implies that 
Fr^^ = 0, and thus the flow of Fr^'^ vanishes. Moreover, the equation for the proper two-point 
function involves, on the right hand side, only the term containing Fr , 

5^^* a\ ^Sfc^^^^^ [drRr]GrU , (96) 

see also Fig.[T21 The term in parentheses stands for the regularised line. We supplement Eq. ((M]) 

with the flow equation for Fr'^\ which, for (f>a = 0, is drawn in Fig. [121 This system of equations is 
still exact. For practical computations it needs to be closed which can be achieved by truncation 
or by supplementing it with equations for one or more higher n- vertices truncated at some higher 
order. Here we truncate by neglecting, in the equation for Ft'^\ the term involving Ft^\ 

drF^Xd = -\jjyrl%^fGrjA%,h} {Gr[drRr]GrW + P{a, b, C, d) , (97) 

where P implies a sum corresponding to all permutations of its arguments. In this way we 
obtain a closed set of integro-differential equations for the proper functions up to fourth order. 
As we will show in the following, they allow to derive, for a particular cutoff time r, a set of 
dynamic equations describing the time evolution of the two- and four-point functions up to 
time t = T. We emphasise that the only approximation here is the neglection of the six-point 
vertex, see Fig. [121 



3.7.3 Dynamic equations 

For the sharp cutoff i?r chosen here the flow equations can be analytically integrated over t. 
As pointed out above, our cutoff implies the connected two-point function to vanish at times 
greater than r, i.e., it can be written as 



Gr.ab — i 



F. 



(2) 



(98) 



where 0{t) evaluates to for r < and to 1 elsewhere, and where ta is the time argument 
corresponding to the field index a, etc. Hence, the precise way in which the cutoff Rr.ab diverges 

a.tta ~ tf, > T is chosen such that /"r +Rt is the inverse of —iGt for all times ta, see Eq. (|95|) . 
Using Eq. ([551) o^^c finds that 



iGr[drRr]Gr)ab = -iGr.abdrW - ta) 0(t - tf,)] 



(99) 

Note that we have not used the specific choice (|92p to derive (|99p but simply the property 
([98ll . i.e., the suppression of any propagation for times t > t. After inserting Eq. ([99l) into 
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Eqs. ([96|l and ([97|l wc can integrate over r and obtain, after some algebra, the integral equations 
determining the flow of the proper functions from Iq to some final time t, 

^ 1 4 

= 9 / ^r,a,acbd'Gr,d.dc, (100) 



^ r,ah 



T.abcd 



'-0 



to,C 
t 



<'L,a6e/G.,,./,r(J^,.,,,,G.^.j.e + (a - c) + (a - d). (101) 

See Ref. [140j for more details of the derivation. Double indices imply sums over field compo- 
nents, spatial integrals and time integrations over the CTP C, from to to t and back to to. We 
furthermore introduced 

Tab = niax{<a, tb}, Tabcd = max{ta, tb, ic, td}- (102) 

The brackets denote terms with the respective indices swapped. 

From Eqs. (|100[) and (|101[) it is clear that for the two- and four-point functions to be defined 
at T = i we need to specify initial functions at t = to- We point out that, within the truncation 
scheme chosen above, we can insert any set of proper two- and four-point functions, as long 
as we set all n- vertices for n = 1,3, and n > 4 to vanish. Hence, the scheme corresponds 
to a Gaussian initial density matrix. Here, we choose the respective classical proper functions 
defined by S in Eq. ([551) . Hence, the initial two- and four-point functions entering Eqs. (jlOOp 
and (jlOip read 

= Sl^^ = iG-^l,, (103) 

^to!abcd = ^abcd = ~{'29/-^)iSabScd + SacSbd + SadSbc)5ciXa ~ Xb)Scixb ~ Xc)SciXc - Xd).{10'i) 

In order to arrive at a set of dynamic differential equations, we finally rewrite Eq. (|100p as 

'^Go.lcGr^t.cb = iScab - 9 / ^T^^.adce ^r^.M GT,t,cb, (105) 
^ Jto.C 

with ScMb = SabScixa - Xb) = SabScita - 4)^(xa - Xf,). Eq. (|105p is the dynamic (Dyson- 

f 2) 1 

Schwinger) equation for the connected two-point function Gr^^^ab- Since = iG^ repre- 

sents a differential operator, the solution of Eq. (|105p finally requires another set of boundary 
conditions to be specified, depending on the form of the differential operator. In our case, this is 
the initial two-point function, fixed by the one-body density matrices as well as Bose statistics, 
see Sects. 13.41 and 13.8.11 We also point out that we have set t = Tcb since the flow of Gr.cb 
stops at the maximum of the time arguments tb,tc. This can be proven from the structure of 
Eqs. poip . (|105p but is more easily seen from the definition Once the hard cutoff r has 
passed the largest time appearing in a (connected) correlation function, the flow stops since the 
forward and backward parts of the CTP over greater times cancel identically in the functional 
integral. 

Let us assume that t = ta denotes the present time, at which Eq. p05[) determines the 
further propagation of Gt^,ab for tb < ta (for tb > ta the solution is then flxed by symmetry). 
All time arguments of the functions occurring on the right-hand sides of Eqs. (jlOOp . (jlOip . and 
therefore all cutoff times Tab and Tf,fgh are evaluated at times t' < t. Hence, in accordance with 
causality, Eqs. (jlOip and (|105p for a given initial Gto,ab{to,to), can be solved iteratively. 

This concludes the derivation of a closed set of dynamic equations for the two-point corre- 
lation function as obtained from a functional RG approach with a cutoff in real time. In the 
remainder of this section we will concentrate on rederiving the dynamic equations obtained 
from the 2PI effective action in NLO of the 1/Af expansion. 



3.7.4 From RG to 2PI next-to- leading order (NLO) 1/7V 

The mean-field approximation, i.e., the HFB equations, cf. Sects. [2731 and [376.31 are obtained 
by neglecting the flow of the four-point according to Eq. (|10ip . F^^^ = F^^K As a consequence, 
the flow parameter of Gr,ab is, for the chosen cutoff, fixed to r = Tab and can thus be neglected. 
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{tf, tg= ?„...min(r,0 ) 




( tg = f„... max(f,„ g ) 

Fig. 13. (color online) The upper equation is the s-channel projection of Eq. (jlOlfl . The second equation 
defines the resummed vertex appearing in the Dyson-Schwinger equation derived from the NLO l/A/" 
approximation of the 2PI effective action [29]. The two definitions are identical in every order of a 
perturbative expansion (see text). Dashed lines denote the s-channel part of the bare vertex F^^^^i,, see 
text. All other symbols correspond to those in Fig. [I2l Letters on internal lines indicate summation 
over field indices and integration over space and time (along the OTP from to to t and back). The 
integration intervals are given in parentheses. 

As a first step beyond mean-field we consider the truncation in which the s-channel scat- 
tering diagram is included beyond the mean-field limit in the loop integral on the right-hand 
side of Eq. (jlOip . The result of this section will be that the obtained equations are equivalent 
to the 2PI equations in NLO of a 1/Af expansion. This truncation corresponds to keeping only 
one channel of each, the classical vertex and the one-loop integral term. The four-vertex then 
only depends on two space-time variables and field indices, ^t^^^bd ~ ^t'ab^<^<^^bd, and enters, 
in Eq. (|105p . as a cutoff-dependent self energy Sr^^^ab = I'r'^^^ab^'^abMb- We can always write 
t = min(t, t'), with t = f' = Tab, for the cutoff parameter in r^'^^^ , see Fig. 1131 upper panel. Since 
the parameter in the (green) vertices on the right hand side is also the minimum of the maxima 
of integration times in the adjacent loop and the respective external time t or t' , one can iterate 
the integral equation in order to obtain a perturbative series of bubble-chain diagrams consist- 
ing only of classical vertices and full, cutoff-dependent propagators. This procedure provides us 
with a proof that St^^ ab is identical to the NLO 2PI 1/J\f self energy shown in the lower panel 

of Fig. m 

A detailed formal proof of this result can be found in Ref. [140] . Here, we provide an 
argument that the above identity can be inferred in a comparatively easy way from the topology 
of the different terms in the flow equations for the two-, four- and six-point functions: Consider 
the untruncated set of equations as displayed in Fig. 1121 First, non- s-channel contributions do 
not generate bubble-chains of the form shown in Fig. 1131 Second, Fr^^ is one-particle irreducible. 
Its contribution to the flow of J"i^' does not give rise to bubble-chains, even if inserted recursively 
into the first diagram on the right-hand side of the flow equation for I^t^^ . In turn, by dropping 
the second diagram with the six-point function and using the s-channel truncation, the iterated 
flow equation generates only bubble-chain diagrams with full propagators as lines. Hence, a 
T-integration of this set of flow equations leads to dynamic equations which include all bubble- 
chain contributions and are therefore equal to those obtained in the NLO 1/Af expansion of 
the 2PI effective action, sec Sect. 13.631 We emphasise that the above topological arguments are 
generally valid when comparing resummation schemes inherent in RG equations of the type of 
Eq. ([94]) . with those obtained from 2PI effective actions. This applies, e.g., to equilibrium flows 
|142I143I144] and thermal flows [1431148] . For a comparison with 2PI results see Ref. [148] , for 
the interrelation of 2PI methods and RG flows Ref. [144j . 
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3.8 Dynamics of a one-dimensional Bose gas 

In this section we apply the theoretical methods summarised above to describe the equilibration 
dynamics of a uniform ultracold gas of bosonic sodium atoms which are confined such that they 
can move in one spatial dimension only, see Sect. 12.4.31 The gas is assumed to be initially far 
from thermal equilibrium. The ensuing equilibration process is found to happen on two different 
time scales. A fast dephasing period leads to a quasistationary state which shows certain near- 
equilibrium characteristics but is still far from being thermal. After this, the system approaches, 
within an at least ten times longer period, the actual equilibrium state. 

3.8.1 Initial conditions 

The 2PI effective action approach is convenient for situations, where at time t = to one has a 
Gaussian state, i.e., a state, for which all but the correlation functions of order one and two 
vanish. In the following we will consider a one-dimensional uniform system, for which the two- 
point functions Fab{x,y) and pab{x,y) are spatially translation invariant. We therefore work in 
momentum space, where the kinetic energy operator is diagonal. Moreover, we choose the field 
expectation value </> to vanish initially. Then, for reasons of number conservation, the equations 
of motion (jTS]) and (|74p will conserve (/) = for all times, see Sect. 12.11 We note that, since 
there is no spontaneous symmetry breaking in one spatial dimension at non-zero temperature, 
the field always approaches zero eventually, irrespective of its initial value. 

Having prescribed initial values i^Qf,(0,0;p) and pab(0, 0;p), the coupled system of integro- 
differential equations (|74p yields the time evolution of the two-point functions, in particular, of 
the momentum distribution 

n{t, p) = ^ (Fh (i, t; p) -f F22 (i, t; P) - a) , (106) 

cf. Eq. ([5]). We will later, in Sect. O distinguish between quantum and classical statistical 
evolution. For the quantum gas, one has a ~ 1 from the Bose commutation relations, while, 
for a gas following classical statistical evolution, a = 0. 

We choose, at f = 0, a Gaussian momentum distribution 

n(0,p) = ^e-P'/"'. (107) 

which constitutes a far-from-equilibrium state if the interactions are non-zero and the corre- 
sponding interaction energy is much larger than the kinetic energy. 

The initial pair correlation function vanishes, [Fii{t, t;p) — F22{t, t;p)]/2 + iFi2{t, t;p) = 0, 
for t = 0, in accordance with total atom number conservation at non-relativistic energies, see 
Sect. O Hence, 

Fn (0, 0; p) = F22 (0, Q-p) = n[Q,p) + a/2, F12 (0, 0; p) = F21 (0, 0; p) = 0. (108) 
As far as the spectral functions are concerned, the Bose commutation relations (j57p imply: 

pil{t,t-p)^ P22{t,t-p) = {), ~ pi2{t,t-p)= p2l(t,t-p) = l. (109) 

Let us investigate the dynamic evolution of a ID Bose gas of sodium atoms with mass 
m = 22.99 u (^^Na) in a box of length L = NgUs, with periodic boundary conditions. We choose 
the numerical grid such that it corresponds to a lattice of Ns points in coordinate space with 
grid constant Cg, and the momenta on the Fourier transformed grid arep„ = (2/as) sm{mT/Ns) 
as explained in detail in Ref. [39]. The results presented in the following are obtained using 
Ng = 64 modes on a spatial grid with grid constant = 1.33 pm. We consider a line density of 
the atoms in the box of ni = 10^ atoms/m. In this case the atoms are weakly interacting with 
each other, such that g = giD = fi^^ni/m, with the dimensionless parameter 7 = 1.5 • 10~^, see 
Sect. 12.4.31 The width of the initial momentum distribution is chosen to be it = 1.3 • 10^ m^^. 
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Fig. 14. (Color online) Momentum-mode dis- 
tribution n{t;p) for the initial state (red filled 
circles, interpolated by red solid line) and 6 
subsequent times t until no change can be ob- 
served for t > 0.6 s. The interpolation of the 
final distribution is shown as a black solid 
curve, n is normalised to the total number of 
atoms, niL — 853. The gas is in a far-from- 
equilibrium state initially, characterised by a 
Gaussian distribution n(0;p), Eq. (|f 07p . with 
width a = 1.3-10^m~^. It is weakly interacting, 
7 = 1.5 -10"^. Since we consider a homogeneous 
gas and a symmetric initial state, the occupa- 
tion numbers are invariant under p —p. 

3.8.2 Equilibration of the quantum gas 

To solve Eqs. ((74|). with the self-energies given by Eqs. ([84|) and the initial conditions in the pre- 
vious section, a parallelised Runge-Kutta solver has been implemented and used on a cluster of 
3 GHz dual processor PCs with up to one node per momentum mode. The correlation functions 
Fabit,t';p) and pab{t,t';p) were propagated, for fixed t', along t, using a second-order Runge- 
Kutta algorithm. After each Runge-Kutta step, the If,p integrals were updated according to 
Eqs. ([85|) . The dynamic equations derived from the 2PI effective action arc, by construction, 
number and energy conserving. While number conservation, by virtue of the 0(2) symmetry of 
each diagram, is given exactly, energy conservation may be violated by the chosen discretiza- 
tion along the time axis. Hence, in order to ensure optimal energy conservation numerically, 
a fourth-order Runge-Kutta algorithm was employed for the propagation of the correlation 
functions along the diagonal t = t'. 

Fig. [14] shows, as a (red) solid curve, the initial Gaussian momentum distribution of the 
gas, on a logarithmic scale, where it forms an inverted parabola. The filled circles indicate the 
numerically calculated modes pi. In the same figure, the time evolution of the distribution is 
shown for different times between t = 0.1ms and 0.6 s. For times greater than about 0.15 s, 
there is only very little change observed. As a function of time, the evolution of the single mode 
occupations is shown in Fig. [151 We observe that the system very quickly, after about 5 /iS, 
evolves to a quasistationary state, and that the subsequent drift to the equilibrium distribution 
takes roughly ten times longer. In passing we note that the mean-field Hartree-Fock (HF) 
approximation, for which 2J^ = XJp = in Eqs. ()74|) . conserves exactly all mode occupations 
and no equilibration is seen. 

In order to estimate to which extent the final distribution approaches that of the actual 
equilibrium state of the gas, we fitted the distribution to the Bose-Einstein-like form n(t]p) = 
[exp{{uj{p) — ji) /kBO{t]p)} — with a p-dependent temperature variable 0{t;p). Here uj{p) 
was derived from the time-derivatives of the statistical function F{t,t';p) at t = t' . If a Bose- 
Einstein distribution is approached the temperature can be obtained from the slope of log(n~^ + 
1) and the chemical potential fj, from its value at cj = 0. Fig. [TBI shows 0{t;p) for t = 0...0.6s. 
Obviously, during the quasistationary drift period, no temperature can be attributed to n{t;p), 
while, for large t, becomes approximately p-independent. We deduce an approximate final 
temperature from 0(0.63; 128/L) = T = 0.35nK with fj, = 1.08 (7lD?^l for the above given values 
of giD and rii, which deviates from the HFB result fi = giD^i by about 8%. 

In summary, the dynamic equations derived from the 2PI effective action expanded beyond 
the mean-field approximation can be applied to derive the evolution of a one-dimensional Bose 
gas starting with a momentum distribution far from equilibrium. Their solutions show a non- 
secular evolution towards a state characterised by a thermal momentum distribution. The 
evolution is characterised by a fast initial dephasing followed by a long slow drift to thermal 
equilibrium. 

A remark is in order: As pointed out in Sect. l2.4!3l a homogeneous ID Bose gas with contact 
interactions represents an integrable system. Hence, at large times, recurrence phenomena are 
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Fig. 15. The normalised momentum-mode oc- 
cupation numbers n{t;p)/niL, corresponding to 
those shown in Fig. 1141 as functions of time. 
Shown are the populations of the modes with 
p = p, = 2Ns/Lam{iTv/Ns), i = 0, 1, iVs/2, 
and one has n{t; —p) = n{t;p). A fast short-time 
dephasing period is followed by a long quasista- 
tionary drift to the final equilibrium distribution. 
Notice the double- logarithmic scale. 



Fig. 16. Momentum and time dependent temper- 
ature variable 0{t; p) obtained by fitting the distri- 
bution n{t;p) = [exp{(ci;(p) — fi) /kB0{t; p)} — 
to the distribution obtained from the results shown 
in Fig. 1151 for different, equally spaced times be- 
tween t = and t = 0.6 s. One observes that, dur- 
ing the quasistationary period, 0.01 s< t < 0.1s, 
no temperature can be associated to the distribu- 
tion. Only at very large times, O becomes approx- 
imately p-independent. 



expected which therefore prevent a thcrmalisation of the system. This does not imply that 
the above results obtained for a finite-time evolution are invalid. We point out that they were 
obtained for a weakly interacting gas for which the approximations made are expected to be 
justified. Nonetheless, questions remain open — whether in an actual system a quasi thermal 
state is reached and, if so, on which time scales and under which conditions the integrability of 
the system drives it away from there again |77ll27lli5] . 



4 From nonequilibrium to kinetic equations 

Far-from-equilibrium dynamics is qualitatively different from evolution near equilibrium where 
distinct properties of the equilibrium state are still "felt" by the system. We have mentioned a 
few times above that the statistical correlation function F is, near equilibrium, no longer inde- 
pendent from the spectral function p. Their interrelation is a manifestation of the fiuctuation- 
dissipation relation well-known in nonequilibrium statistical physics. Before we set our focus 
on how near-equilibrium behaviour emerges during the dynamical evolution starting far from 
equilibrium we illustrate the fluctuation-dissipation theorem in the framework of linear-response 
theory. 



4.1 Linear-response theory 

Analysing a simple quantum mechanical linear dissipative model, such a theorem can be easily 
derived from the properties of perturbation and response functions. Let us consider a general 
many-body Hamiltonian H, perturbed in time around some equilibrium Hamiltonian Hq by 

Hiit) = -/ie-'('^+")* + h.c, (110) 

where h is an operator inducing a small perturbation, and the regulator suppresses the pertur- 
bation for t ^ — oo. In linear- response theory one studies the fluctuation 
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where X/t./i is commonly called the dynamic polarisability. In perturbation theory, one finds 

{m\f^n){n\h\m) (m|/i|7i)(?i|/t|?7i 



(112) 



with LOrnn = ^ ^o|?i) = £'Ti|?^), Z = cxp{— /3i?„}, and where thermal equilibrium is 
assumed for t ^ — oo. One also defines the dynamic structure factor 

Sficu) = ^e-''^"|(n|/|m)|2j(c. - c.„„), (113) 

m,n 

which vanishes for lu < and T ^ 0. In the case that f ~ h one then finds that the dynamic 
structure factor is related as follows to the response function Xf = X/t./- 



X/(w) = - / duj' 



OO 



uj — uj' + ie LU + Lu' + ie 



(114) 



Using the relation (x — a ± ie)~^ = ■p(x — a)^^ T ii^6{x — a) involving the Cauchy principal 
part V one finds that the imaginary part of the response function is related to the dynamic 
structure factor as follows, 

lxnxf{oj)=Ti{l-e-PnSf{i^). (115) 

where we have used that, at finite temperatures, one has Sf{uj) = e'^'^Spi-uj). The imaginary 
part of X is usually called the spectral density related to the operator /, 



Pf 



(uu) = ^-1^(6-^*^'" -e-'='^")|(n|/|m)|227r(5(c.~L.„„0 = 2Imx/(^). (116) 



In second-order perturbation theory one finds the rate of energy transfer to the system to be 
dE/dt = 2a;Imx/(w) + terms oc exp{2ia;} which provides us with a physical interpretation of 
this quantity. 

A sum-rule argument can now be used to relate the zeroth-order moment of the structure 
factor ()113p to the fluctuations of /, 

dLo{Sfiu) + Spico))^{{f,P})p, (117) 

where = Z"^ e~^^'^{m\ ■ |to), and where the completeness of the states \m) has been 
used. Using again Sf^cu) = e^^ Sj:i{^Lj) one finds the fluctuation-dissipation relation 

({/, '^^Pfi^) ("be(c.) + ^) , (118) 

with the Bose-Einstein distribution ^,-^-^{10) = (exp{/3ct;} — 1)~^. 



4.2 Equilibrium field theory 

In thermal equilibrium, correlation functions can be determined using the functional quantum 
field theoretical approach discussed in Sect. [3] if the real-time closed time path is replaced by 
the imaginary path C ~ [0, (3 denotes the inverse temperature. Let us consider a spatially 
uniform system. Since thermal equilibrium is time translation invariant, the two-point functions 
depend only on relative space and time coordinates. In Fourier space one has, e.g.. 
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Fig. 17. (Color online) Left panel: Ratio of the envelopes of the unequal-time correlation functions, 
^{t-p) = [{Fii{t,0;pf + Fi2it,0;pf)/{puit,0;pf + pi2it,0;pf)]'-^^/n{t;p), for four different momen- 
tum modes, as a function of time. The evolution corresponds to that shown in Fig. 1151 Sect. 13.8.21 
Due to the normalization with respect to n{t;p) all ^{t;p) are of the same order of magnitude. ^ is 
a measure of the interdependence of the statistical and spectral functions, and its settling to a con- 
stant value during the quasistationary drift period indicates that these functions become, as in thermal 
equilibrium, connected through a fluctuation-dissipation relation. Right panel: Evolution of the kinetic 
and interaction contributions to the total energy of the gas. During the quasistationary drift these 
contributions assume the same order of magnitude, calling in mind the virial theorem. 



Since, in equilibrium, the path integral determines the expectation value with respect to the 
density matrix p cx exp{— the initial and final states must be set equal. This leads to the 
periodicity condition for tlie imaginary-time propagator, Gab{x,y)\xQ=o = Gabix,y)\xo=-ip- 
Inserting Eq. ([5]) into this condition one obtains the relation 

Fir^(.;,p) = ^ (nM^) + pif (^^P) (120) 

which, integrated over uj, is identical to the fluctuation-dissipation relation (jllSp |20ll50j . In 
deriving Eq. (jl20p we have used that the time xq = is the earliest on the imaginary path C 
while xq = —i(3 is the latest, such that sgn^(a;o — t/o) yields opposite signs in the two limits. An 
equivalent fluctuation-dissipation relation can be derived, using the same arguments, for the 
equilibrium self-energy contributions Sp and Sp. 

We emphasise that in nonequilibrium fleld theory, the relation (|120[) no longer holds. Far- 
from-equilibrium dynamics in this respect allows for a greater variety of solutions. Nonetheless, 
nonequilibrium fleld theory encompasses the near-equilibrium as well as vacuum theories, where 
"-be(w) 0. 



4.3 Near-equilibrium time evolution of a ID Bose gas 

We now study the equilibration process of the one-dimensional Bose gas presented in Sect. 13.8^ 
with respect to the characteristics of near-equilibrium evolution introduced in the previous 
section. For this we consider the time dependence of the ratio of the envelopes of the unequal- 
time correlation functions, speciflcally, £,{t;p) = [(Fii(i, 0;p)^ + Fi2(t, 0;p)^)/(/9ii(i, 0;p)^ + 
Pi2{t,0;p)'^)]^^^ /n{t;p). Confer Refs. |46|48j for more details. Similar studies for relativistic 
gases have been presented in Refs. [1511152] . Fig. [T7] (left panel) shows ^, for four different 
momentum modes, as a function of time. Due to the normalization with respect to n(t]p) all 
^{t;p) are of the same order of magnitude. However, they show a distinct time evolution during 
the dephasing period, before they settle to a constant value during the quasistationary drift. ^ 
is a measure of the interdependence of the statistical and spectral functions, which, in thermal 
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equilibrium, arc connected through the fluctuation-dissipation relation given, in the momentum- 
frequency domain, in Eq. ()120|) above. Hence, the time at which ^ becomes stationary indicates 
that F and p are linked to each other long before the momentum distribution becomes thermal, 
as can be seen from Fig. [TBI 

A further signature is found when comparing the kinetic and interaction contributions to 
the total energy as shown in Fig. 1171 right panel. During the quasistationary drift, these con- 
tributions are constant and of the same order of magnitude, calling in mind the virial theorem. 

In summary, during the drift period, as can be seen in Figs. [Ml and I16[ the system is not 
yet in equilibrium as far as the momentum distribution and temperature are concerned. The 
statistical and spectral correlation functions nevertheless are already locked to each other as 
shown in Fig. I17( left panel. Hence, one may expect that during this period the dynamics can 
be described using kinetic or transport equations, i.e., that it fulfills the requirements for the 
approximations implied by quantum Boltzmann equations. We shall study this in more detail 
in the following sections. 



4.4 Transport equations 

So far we have dealt with initial-value problems that describe the time evolution by means of 
in general coupled equations of motion for time dependent correlation functions, given specific 
values for them at the initial time. Alternatively, the time-evolution of many-body systems 
is very often described in terms of kinetic or transport equations [20]. In general, the aim of 
quantum kinetic theory is to find evolution equations for distribution functions /(x, t) or /(p, t), 
interpreted, e.g., as particle number densities in x- or p-space. These distribution functions can 
then be used to derive various transport properties of the many-particle system as, e.g., current 
of charge or energy. For that reason, the evolution equations for / are generally referred to as 
transport equations. 

Kinetic descriptions usually neglect the effect of correlations between different times of the 
evolution, i.e., they build, to a certain extent, on a Markovian approximation. In particular, 
they neglect the initial dynamics directly after a change in the boundary conditions which drive 
the system out of equilibrium. This shortcoming is cured in a dynamical approach as discussed 
above. The buildup of correlations beyond the kinetic approximation, in these equations, is taken 
into account by means of non-Markovian integrations over the evolution history of correlation 
functions. 

In the context of non-relativistic systems, in particular cold atomic gases cf., e.g., Refs. pp] . 
jl58ll54ll55ll56ll57ll58ll59ll60l45l44) . See also Refs. [lf)lllf)2ll68ll64lie>5ll()f)ll()7ll68ll6Hj . 
|170|171|172j in the context of relativistic physics. A comparison of dynamic equations with their 
kinetic approximation for relativistic dynamics has been given in Refs. |161I173I174I175I171I172| . 
In this section, we derive transport equations from the 2PI dynamic equations (|74p . in leading 
order (LO) and next-to-leading order (NLO) of a gradient expansion. See Refs. [172|49j for 
more detailed discussions. The field expectation value is set to zero as before. If not stated 
explicitcly otherwise, field indices a,b,... are suppressed in the following, and all products of 
correlation functions are to be taken as matrix products. 

We rewrite the exact dynamic equations ([7^ with the self energy given in NLO 1/Af ap- 
proximation as described in Sect. l3.6!3l in order to get evolution equations with respect to centre 
and relative coordinates, providing a starting point for a gradient expansion and transforma- 
tion to Wigner space. For simplicity we consider a system with spatially homogeneous initial 
conditions for the two-point function. The initial values of the spectral function are fixed by 
the commutation relations, cf. Eq. (|109p . /9o.a()(x, y) = --ia2.abS{x. — y). 

To meet the requirements of the Fourier transformed equations in Wigner space on the one 
hand and to benefit from the spatial homogeneity of the system on the other hand, we rewrite 
the equations in terms of relative and centre space-time coordinates 

X i^i^,s = x-y ^ X = X + s/2, y = X - s/2 (121) 
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where the two-point functions are to be transformed as, e.g., F{x,y) i— > F(X,s), Sp{x,z) t-^ 
EP{X + s'/2,s- s'), F{z,y) ^ F{X + (s' - s)/2,s')- We additionally introduce s' = z - y 
which later on serves as the integration variable in the memory integrals. 

Using these definitions in Eqs. ((74)) . and adding these equations to the corresponding equa- 
tions for the time derivatives with respect to yoi one obtains differential equations with respect 
to the centre time coordinate X^: 



dxoF{X,s) = iaJ \Hib{X + - H^X - + S^^\X + - S^"Hx '-)]FiX,s) 



+ / e{Xo + 4 



So. 



SR{X + -,s-s')F{X 



G''{X+'-,s-s')S^{X 
F{X+'-,s~s')S^{X^ 



t^,s') + S^iX+'-,s- 



s')G'^{X 



(122) 



dxoPiX, s) ^ ^<J2{[HMX + £) _ Hib(X - ^) + S("HX + - S^'Hx ^)]p{X, s) 



[S^X+'^,s 



s')p{X 



p{X + -,s-s')S^{X 



.,s') + SP{X + -,s- 



s')G^{X 



G"{X 



s')SP{X 



^')]}- 



2 ' ' 2 

Note that, introducing retarded and advanced Greens functions and self-energies, 

G'^ix, y) = 0{xo - yo)p{x, y), G^(x, y) = -0{yo - xo)p{x, y), 
y) = d{xo - yo)2J''{x, y), S^{x, y) = -e{yo - xo)S''{x, y), 



(123) 



(124) 
(125) 



allowed us to send the integration limits of Sg ii^ Eqs. (|122p . (jl23p to ztoo. While the above 
equations have been obtained by adding the equations for dxoF(x,y) and dygF{x,y), etc., a 
second set of equations for the derivatives of F{X, s) and p{X, s) with respect to sq results when 
subtracting the respective expressions. These equations are provided and discussed further in 
Appendix A of Ref . [49] . 



4.4.1 Approximations 



With the aim to derive transport equations, we apply the following approximations: 



(j) The ^-function is neglected in the evolution equations for F, Eqs. (|122[) . taking into 
account that the correlations disappear for large relative times. This corresponds to sending 
the initial time tg to the infinite past. Note that, since an interacting system could have reached 
equilibrium at any finite time, transport equations are initialised by specifying F and p at a 
finite time using the equations with to — > — oo as approximate description |172j . 

(ii) A gradient expansion is applied with respect to the centre coordinates X. In leading 
order this expansion corresponds to a Markovian approximation which neglects the variation of 
the correlation functions along the centre time coordinate. The approximation is valid as long 
as this variation is small on the time scale set by the size of the memory kernel. Physically, non- 
Markovian effects are expected to be small for systems not too far from thermal equilibrium. In 
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higher orders, this expansion is a derivative expansion of the correlation functions with respect 
to the centre time and takes into account non-Markovian effects. 

In leading order of the gradient expansion, one obtains from Eqs. (|122p and (|123p . by keeping 
the dependence of all functions S and G on X only and combining retarded and advanced 
functions to spectral functions, the evolution equations 



(126) 



dxoFiX, s) = ia2 J [f{X, s - s')SP{X, s') - p{X, s - s')S''{X, s') 
dxoP{X,s) = 0. (127) 

Note that we have assumed homogeneous initial conditions such that all functions remain 
constant along the spatial centre coordinate. In next-to-leading order in the gradient expansion 
with respect to the centre time coordinate, these equations receive additional corrections. They 
are given explicitely in Ref. [49) , and we suppress them here for conciseness of the presentation. 
The finally resulting NLO equations will be quoted below. 



4.4.2 Transformation to Wigner space 



The transformation to Wigner space involves a Fourier transformation with respect to the 
relative coordinate s 

F{X,p) = J^e^P^FiX,s), E''{X,p)= j^e^P^E\X,sl (128) 

etc., where ps = pqSq — p • s. 

For a spatially homogeneous system one finds that the diagonal matrix elements of F are 
purely real, while the off-diagonal matrix elements are purely imaginary. 

We apply the transformation to the previously derived equations of motion and obtain 

Vx„[F{X,p)] = ia2{FiX,p)SPiX,p) ~ p{X,p)S^iX,p)} (129) 
VxMX,p)] = (130) 

as the leading-order transport equations. In deriving these one needs to take care of the integra- 
tion limits when interchanging the derivative with respect to the centre time with the relative 
time integration. This enters the definition of the derivative operator 

/2Xo r 
dso / d^se'P'dxoFiX,s) 
-2Xo J 

= dxoF{X,p) - 2 [e^'PoXopf^x, so - 2Xo, p) + e-^'P'>'''' F{X, so = -2Xo, p)] . (131) 

Assuming lims^^oo F{so) = = limsQ_>oo p(so) yields limxo-*oo I^Xo ~ 9xo- The contributions 
from the integration limits account for the fact that the correlation functions are initialised at 
some time Xq = yo = 0. These contributions can be removed for sufficiently late times as the 
time correlations are expected to vanish for sufficiently large relative times. 
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In ncxt-to-lcading order of the gradient expansion, the transport equations receive additional 
contributions as follows. 



Vx, [F] = LO + ^2 1 [dxoS'^°'^ [X)] dp,F{X,p) 

-l\{E+iX,p), F{X,p)}o + {S^{X,p),p+{X,p)}o 



VxAp] = <J2^[dxoS^°\x)]dpjix,p) 

-^[{S+{X,p), p{X,p)}o + {S'iX,p), p+{X,p)}o 

These contributions involve Poisson brackets with respect to Xq and pq: 

OA dB dA dB 

{A,B}o = 



(132) 



dpn dXo dXo dpo ' 



(133) 



(134) 



4.5 Comparison of dynamical and transport equations 



It is desirable to compare the transport equations derived in the preceding section, in leading 
order, Eq. (|129p . and next-to-leading order, Eq. (|132p . of the gradient expansion, with the 
equations of motion (|74p . This is achieved in two steps. First, one calculates F and p using 
equations ([74]) in the NLO l/A/'-approximation. The results of this calculation are taken to 
calculate, after a Fourier transformation to Wigner space, both the left-hand and the right- 
hand sides of equations (|129p and (|132p . The left-hand side yields the time derivatives of the 
solutions of the 2PI dynamic equations with respect to the centre time coordinate, while the 
right-hand side generates the corresponding approximative derivatives in leading order and 
next-to-leading order of the gradient expansion. 



4.5.1 Physical setup 

We consider the same dilute homogeneous one-dimensional (D = 1) gas of sodium atoms as 
defined in Sect. 13.8.11 For the results presented below the gas was chosen to be confined in a 
periodic box of length L = NsQs ~ 43 pm. The spacing of the numerical grid is = 1.33 pm, 
and the number of lattice points is Ns = 32. Using different initial line densities ni between 
10^ m~^ and 10^ m~^ the total particle number varies in the range from 4 to 400 particles. 

Three different interaction strengths 7 are considered, see Sect. 12.4.^ 7 = 1.5-10^'^, 0.15, 15, 
with = mgni kept fixed by choosing the line densities accordingly. As will be discussed in 
Sect. [51 the effect of quantum fluctuations is, in this way, increased with 7 while the dynamics in 
the classical statistical approximation remains unchanged, see Sect. [5] The correlation functions 
are initialised as described in Sect. I3.8TT1 i.e., the initial momentum distribution i^(0, 0,p) is 
chosen to be Gaussian, peaked around p = 0, with a width of cr = 6.5 • 10^ m~^, and the Bose 
commutation relations fix the initial spectral function, p(a;o,a;o,p) — —ia2- 



4.5.2 Numerical solution of the dynamic equations 

The 2PI dynamic equations are solved using the techniques described in Sect. 13.8.21 and their 
solutions transformed to Wigner space. We expect the transport equations to reproduce, to 
a good approximation, the full dynamical evolution, if the two-point functions vary slowly 
with respect to the centre time coordinate Xq as compared to the change with the relative time 
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Fig. 18. The statistical two-point function F for two different momentum modes as a function of the 
relative and centre time coordinates so and Xo, respectively. Left panel: momentum mode rip = 0, right 
panel; Up = 8. The figures illustrate that the correlation function vanishes for sufficiently large relative 
times. See the main text for the parameters chosen. 
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Fig. 19. (color online) Spectral function 
/5 (green solid lines), and statistical func- 
tion F, normalised to the occupation 
number of the momentum mode p (blue 
dashed lines), as functions of po, for dif- 
ferent times Xo and different momen- 
tum modes Up. The difference between 
the respective spectral and normalised 
statistical functions decreases, indicat- 
ing the emergence of a fluctuation- 
dissipation relation. The good corre- 
spondence for the zero mode rip = at 
the initial time is related to the choice 
of the initial condition. 



coordinate sq. Figure [TSl shows F{Xq, sq, p) for two different momentum modes. We see that for 
sufficiently late times the functions fall off to zero with increasing relative time Sq ■ The spectral 
function p shows a similar behaviour [49j . and the self energy contributions E^'P decay even 
faster in sq as they involve F and p at least to the third power. Compared to their oscillations 
along the relative time direction, only a weak dependence on the centre time coordinate is 
found. The form of the correlators in the temporal plane shown in Fig. [18] is to be compared 
with that of the correlators for an ideal gas which show an undamped oscillatory dependence 
in the sq direction, with the frequency given by the free dispersion po = p^/2to. 

As discussed in Sect. 14. 2] the two-point functions F and p near equilibrium arc expected to be 
related by the fluctuation dissipation relation (I120p if the system is close to equilibrium [201150] . 
As can be seen in Figure [T9] the function p2i{Xo,p) approaches Fii{Xo,p)/ {n{Xo, p) + ^) during 
the time evolution described by the 2PI dynamic equations ([74]) . where n{Xo, p) ~ Fu {Xq, sq = 

o,p)-i 



4.5.3 2PI dynamical versus transport equations 



We compare the time evolution as derived from the 2PI dynamic equations to their kinetic ap- 
proximations following the procedure summarised in the introduction to Sect. 14.51 The absolute 
values of the left-hand side (LHS), the right-hand side in leading order (LO) and the right-hand 
side in next-to-leading order (NLO) of Eqns. ()129[1 . ()132[) are drawn in Figure [20] 
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Fig. 20. (Color online) Time-derivative of the statistical two-point function as a function of centre 
time Xo for two different momentum modes, Tip = 4 (left panel) and Up = 8 (right panel): Comparison 
of results from 2Pf dynamic equations and those from their kinetic approximation. The (red) solid 
lines correspond to the left-hand side (LHS), the (blue) dashed(-dotted) lines to the the right-hand 
side in leading order (LO) (next-to-leading order, NLO) of the gradient expansion, Eqns. (|129|) and 
(|132[) . respectively. In the left column of each panel, the frequency po is that of the peak of the spectral 
function, cf. Fig. [19] ('on-shell'), in the respective right columns, po has been chosen two half- widths 
away from it ('off-shell'). From top to bottom, the line density ni is rescaled as well as the interaction 
parameter 7, so that mgni = fni is kept fixed: Top row: case A ni = lO^m^^, 7 = 1.5 ■ 10^'^; middle 
row: case B ni — 10^m~^, 7 = 0.15; bottom row: case C ni = lO^m^^, 7 — 15. 



As already seen in Fig. [TS] we observe three generic time regimes |29| : Strong oscillations 
characteristic for early times, slow drifting at intermediate times, and a late-time approach to 
equilibrium characterised by vanishing time derivatives. 

Close to the initial time, the strong oscillations of F{Xo, so,p) and p{Xo, so,p) in the centre 
time Xq are due to the finite integration limits sq = ±2Xo in the sq direction. The intermediate 
drifting regime is reached (for 7 = 1.5 • 10"'^ at Xq ~ 0.003 s) when the contributions of the 
integration limits can be neglected. Non-vanishing derivatives of F and p with respect to Xq 
now solely result from the evolution of the correlations -F(Xo,so,p) and p(Xo,So,p) in Xq. 
As can be seen in Figure [TBI the change with Xq is slow compared to that with sq. Finally, 
equilibrium is approached when all time derivatives vanish. 

We consider in more detail the cases corresponding to (A) weak, (B) moderate, and (C) 
strong effective interactions between the atoms: 

Case A corresponds to weak interactions, 7 = 1.5- 10^^. The evolution is expected to be well 
described by Boltzmann-type equations after the initial oscillations have damped out. Figure 
\20\k (top row) shows the left-hand side (LHS, solid line) as well as the right-hand side of the 
transport equations in leading order (LO, dashed line) and next-to- leading order (NLO, dashcd- 
dotted line) of the gradient expansion for two different momentum modes rip = 4 (left panel) 
and Up = 8 (right panel). In the respective left columns, the frequency po is that of the peak 
of the spectral function, cf. Fig. [12] ('on-shell'), in the right columns, po has been chosen two 
half-widths away from it ('off-shell'). One finds that the 2PI dynamic and the kinetic equations 
in general give the same results once the occupation numbers remain constant. 

The NLO contributions depend on the pq derivatives of the correlation functions. The cor- 
relation functions in Wigncr representation are non-zero for a continuum of frequencies pq, see 
Fig. 1191 This results in non- vanishing NLO contributions. However, as shown by the red dashed 
line, no significant contributions to the transport equations are found after the decay of the 
initial oscillations. 

Case B. We increase the interaction parameter to 7 = 0.15 while the initial line density 
is decreased to ni — 10^. In this way g cx 7711 increases by a factor of 10 such that quantum 
statistical correlations grow in importance, see Sect. [5] Fig. 1201 shows that equilibrium is reached 
faster, in particular for the higher momentum modes, while the intermediate drifting regime 
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Fig. 21. (Color online) Time dependent 
decay constant F as defined in Eq. (|136|l , 
where the time derivative dF/dXo is re- 
placed by the left-hand side (LHS, solid 
line) , the right-hand side in leading order 
(LO, dotted line) and in next-to-leading 
order (NLO, dashed line) of the gradient 
expansion, Eqs. (|129|l and (43) in gS], 
for the on-shell and the off-shell case as 
a function of the centre time. From top 
to bottom, the line density ni is rescaled 
as well as the interaction parameter 7 
with 7711 = const., compare Fig. 1201 The 
three rows correspond to cases A-C as 
described in the text. Colors indicate the 
momentum mode: rip = 4 (red), rip = 8 
(blue). Up = 12 (green). Non-constant 
(and negative) values for F result from 
the fact that the time regime of expo- 
nential approach to equilibrium is not yet 
reached. 



observed in case A is reduced. For the lower momentum modes we get similar results as in case 
A while for the higher momentum modes there is an essential difTcrence to the preceding case: 
correspondence between 2PI and transport equations is not reached until equilibration occurs. 

Case C. We finally choose strong coupling, 7 = 15, and decrease the line density to 
ni = 10^m~^. Qualitatively similar results as in the preceding cases and a faster approach 
to an equilibrium configuration are found, see third line of Fig. 1201 

To study in more detail the late-time behaviour, we assume that, at late times, the statistical 
correlation function decays exponentially to its equilibrium value with a decay constant r{p), 

F{Xo,p) « FiXo = oo,p) + ziF(p)e-^(P)^«, (135) 

where AF{p) is some constant independent of Xq- In Fig. [21] we plot 

d^F{Xo,p)/dX§ 
r{Xo,p) = — ...^ , (136) 
dF{XQ,p)/dXQ 

where OF/OXq is given by the LHS, as well as by the LO and NLO expressions on the RHS 
of Eq. (|129p . We focus the range of times Xq to those where r{Xo,p) is settling to a constant, 
indicating the emergence of an evolution according to kinetic theory. The top row of graphs 
in Fig. [21] shows that the different mode evolutions are settling to an exponential decay at 
times between 0.07 and 0.09 seconds, deep in the drifting regime. Hence, although a fluctuation 
dissipation relation is established almost an order of magnitude in time earlier, the kinetic 
approximation becomes strictly valid only at very late times. Fig. [21] also shows that it is in 
general not sufficient to take into account only the LO approximation in the gradient expansion. 

Let us furthermore study the dependence of the decay constants F on the line density ni. 
Fig. 1221 shows F for five different densities ni. We find an approximately linear dependence 
of F on ni which indicates that the source of damping is rather an off-shell two-body than a 
three-body scattering effect. 
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Fig. 22. (Color online) The decay con- 
stant r as defined in Eq. (|f 36|l as a 
function of the particle density ni, ex- 
tracted from the solution of the 2PI dy- 
namic equations at times tkin as indi- 
cated with (blue) diamonds in the in- 
set figure, after which the evolution can 
be described to a good approximation 
by an exponential decay. In the inset, 
the corresponding times tdrift at which 
the slow drift of the occupation num- 
ber sets in are indicated with (green) 
circles. The error bars indicate the vari- 
ation of r, ikin, and fdrift over the dif- 
ferent momentum modes. The red line 
is a linear fit. 

4.6 Summary 

In this section we have studied further the dynamics of our equihbrating one-dimensional Bose 
gas, thereby focusing on the comparison between a fully dynamical approach on the basis of 
the 2PI effective action and the corresponding kinetic approximation in the form of Boltzmann- 
type transport equations. Transport equations are derived from the 2PI dynamic equations for 
the two-point correlation functions by means of a gradient expansion with respect to centre 
coordinates and a subsequent Fourier transformation with respect to relative coordinates, i.e., 
a Wigner transformation. Furthermore, the details of the initial state are neglected by sending 
the initial time to minus infinity. We found that the time evolution of the correlation functions 
with respect to centre time occurs relatively slowly as compared to the oscillations with respect 
to relative time. 

A comparison of the statistical and spectral correlation functions showed that they begin 
to be locked to each other as predicted by the equilibrium fluctuation-dissipation relation. This 
occurs approximately at the same time as when the transport equations set in to be a good 
approximation to the dynamical evolution according to the full dynamic equations. For weak 
couplings one observes good correspondence of dynamic and transport equations after an initial 
period of oscillations. However, off-shell effects are not covered by transport equations before 
equilibration occurs. Increasing the dimensionless interaction strength parameter 7 leads to 
significant differences between the 2PI dynamic and the transport equations. Hence, the kinetic 
description of the time evolution in terms of an exponential decay with decay rate F generically 
sets in to be valid only at very large times, when no essential change in the momentum profile of 
the system occurs any more. The late-time evolution is predominantly due to two-body off-shell 
scattering effects. 



5 Quantum vs. classical statistical dynamics 

Many experiments with Bosc-Einstcin condensates of dilute gases have shown that the Gross- 
Pitacvskii theory represents, in practice, a very good approximation to describe static and 
dynamic features of these systems. Despite the fact that the first-order coherence reflected by 
this equation has its origin in the quantum nature of the Bose condensation phenomenon, the 
GPE arises as the classical field-theory approximation of the underlying quantum many-body 
problem. It thus neglects all quantum statistical fluctuations contributing to the dynamics of 
the scalar fleld. However, it is the role of these quantum statistical fluctuations which is of 
central importance for our quantitative understanding of quantum many-body dynamics, so far 
the least explored in the dynamical world. 
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Two cases should be distinguished in this context: If the real-time dynamics of a Bose gas is 
dominated by classical statistical fluctuations then it can be well approximated by a large num- 
ber of numerical integrations of the classical field equation (GPE) and Monte Carlo sampling 
techniques. Such simulations are done in many different areas in physics and to some extent form 
a research field in its own. In the context of ultracold Bose gases, see e.g., Refs. |176ll77ll78j . 
Comparisons between simulations of classical field dynamics and the nonperturbative 2PI ap- 
proach described in Scct.[3]have been presented in Refs. |179|37j . Simulations take into account 
nonperturbative dynamics, however, they neglect all quantum corrections. For fermions, a cor- 
responding classical statistical description does not exist. The other case concerns dynamics 
where quantum fluctuations are relevant. In the following we quantitatively determine the role 
of quantum fluctuations for a time-evolving Bose gas. 

Many experiments concerning ultracold Bose gases fall short of being sensitive to quantum 
statistical fluctuations and can be accurately described by classical field theory. This encom- 
passes mean-field approximations which are applicable if classical fluctuations are small as 
well as the case of strong classical fluctuations which requires simulations of classical field 
equations. The importance of classical statistical fluctuations can rise if the gas is sufficiently 
dense. A combination of low densities and strong self-interaction can lead to enhanced quan- 
tum fluctuations as compared to classical statistical fluctuations. As discussed in Sect. 12.41 
zero-energy scattering resonances, particularly the so called magnetic Feshbach resonances, as 
well as lower-dimensional gases so far have played a leading role in the creation of strong in- 
teractions in degenerate atomic quantum gases. Present-day experimental techniques allow for 
resonance-enhanced scattering lengths larger than the mean interatomic distance n~^/'^, such 
that the system is no longer in the coUisionless regime. 

In the following we discuss the difference between the quantum and classical statistical 
theory which can be expressed in terms of interaction vertex terms for the quantum theory which 
are absent in the classical statistical theory. This implies that the classical generating functional 
is characterised by a reparametrisation property, which allows one to scale out the dependence of 
the dynamics on the scattering length a. As a consequence, for the classical dynamics the effects 
of a larger self-interaction can always be compensated by a smaller density. It is shown that 
quantum corrections violate this invariance property. They become of increasing importance 
with growing scattering length or reduced density. On this basis one derives a condition which 
experimenters may use to find signatures of quantum fiuctuations when preparing and probing 
the dynamics of ultracold gases. To illustrate the differences we compare quantum and classical 
evolution for the example of a one-dimensional Bose gas equilibrating in one spatial dimension 
as it has already been the subject of Sects. [3^ and 14.51 



5.1 Functional-integral approach to classical statistical dynamics 

The use of functional methods to describe the dynamics of classical correlations dates back to 
the work of Hopf in the context of statistical hydrodynamics |180j . A field theory for the descrip- 
tion of classical fluctuations in terms of noncommutative classical fields was first suggested by 
Martin, Siggia, and Rose (MSR) [181j and has been extensively used in critical dynamics near 
equilibrium |182| . This theory has been reformulated later in terms of Lagrangian field theory 
employing functional methods [18311841185118611871188] . In these field theoretical approaches 
to classical statistics, a doubling of the degrees of freedom occurs. For example, in the generat- 
ing functional for Green functions, besides each field appearing in the fundamental Lagrangian, 
a second 'response' field is integrated over. The functional integral approach to quantum field 
dynamics developed by Schwinger and Keldysh employs the closed time path (CTP) contour 
[1891190] in the time-ordered exponential integral as introduced in Sect. 13.31 The doubling of 
fields in the MSR and Lagrangian approaches to classical dynamics corresponds to the fields 
evaluated separately on the two branches of the Schwingcr-Keldysh CTP |191|192|193| . Implica- 
tions of the differences between the classical and quantum vertices, similar to the case considered 
in this article, have been discussed, for other theories, e.g. in Refs. |194|195|196|150|197|193|198] . 
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5.2 A quantum-mechanical example 

Let us begin with a simple quantum mechanical example, that of a quantum harmonic oscillator. 
As is well known and was mentioned in our introduction to mean-field theory, the coherent states 
([1]) formed from the energy eigenstates of the oscillator describe classical motion as it is, e.g., 
seen in the trajectory of the mean value of the position operator. The wave packet shows a 
dispersionless motion in the potential. As a reason for this special behaviour one can regard 
the absence of genuinely quantum effects like tunneling or quantum reflection. 

In the context of field theory it is common to regard mean-field approximations, including 
the HFB theory as classical approximations in the statistical sense. The reason for this is 
that the Hamiltonian is at most quadratic in the field operators and their derivatives. Such a 
Hamiltonian can always be diagonalised and represents an essentially free system. 

Before discussing these matters further let us shed a bit more light on the classicality of the 
quantum harmonic oscillator using the path-integral formulation of initial-value problems. As 
discussed in detail in Sect. 13.31 a path integral allowing to derive the time-evolution of operator 
expectation values from a given initial state involves a time-integration over a closed time path, 
leading from the initial time to to the maximum time of interest and back to to. Hence, the 
argument of the exponential in the measure, cf. Eq. (j69p . can be written as 



Sc[v^] = S[^+]~S[^-], (137) 

where (p^ are the fluctuation fields evaluated on the outward (-I-) and backward (— ) branches of 
the CTP, respectively. Note that the relative minus sign stems from the reverted time integration 
on the '— ' path. The CTP provides a doubling of fluctuation coordinates ip as compared to the 
path integral for the simple quantum mechanical transition matrix element (|42l) . In the case of 
the harmonic oscillator with frequency uj the difference (|137[) can be cxplicitely written as 

= \ [^t^o - ^0^0 - Jdt y+ {d? + Lu^) ^+ - if- {d^ + uj^) I , (138) 

where ip^ = dt^p^, and an integration by parts produced the boundary terms involving the 
coordinate lPq = (p^(to) and velocity ip^ evaluated at i = io- Expressing now the coordinates 
evaluated on the -I- and — branches of the CTP in terms of a "centre" coordinate (p and a 
"relative" coordinate (p, as ip^ = ip ± ip/2, we can rewrite Sc as 

Sc[^]^S[ip,^ = -l^ipodtipo+ J dtlp{d^+co^)^Y (139) 

where we have performed two further partial integrations in order to get rid of any terms 
involving derivatives of (p in the integrand. In this way, also the boundary terms combine to 
the single term left in Eq. p39p . Inserting p39p into the path integral ((69|) for initial- value 
problems one obtains 

V^+V^- (^+|po|¥'o")e*^^''"^l"'''^"" = / VipV(pp[ipo,^o]e''^'^-^, (140) 

with p[ipo, ^Po] = {'P'o + <^o/2|po|'/'o ~ '?o/2), since the Jacobian of the transformation to ip and tp 
is equal to one. This path integral is of the form discussed in the context of classical statistical 
physics and contains, besides the integration over the physical coordinate ip, an additional 
integration over an auxiliary coordinate ip, see, e.g., Rcfs. .183 184 185 186 187|188j . The final 
step is done by observing that the integral over Ip can be analytically performed, which gives, 
as integrand of the remaining path integral, a functional delta distribution, 

VfVipp[ipo,f>o]e'^^'''^^ = [ d(^od7roTybo,7ro] / V" ip S[{d^ + lu^) ip]. (141) 
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Here, tt = dt^p denotes the momentum canonically conjugate to the harmonic-oscillator coor- 
dinate Lp, and the initial phase space distribution VF[</joj ""o] is defined as the Fourier transform 
of the density matrix in (/Jo-'?n representation, 

W^o, TTo] = j d<^o Pbo, ^o] exp{-i<^o7i"o}- (142) 

Note that the phase factor in Eq. p42p originates from the boundary terms obtained in the 
partial integrations in Eqs. (|138[) and (|139[) . The double-primed measure Vip excludes the 
integrations over the initial-time coordinate ip^ as well as over the next-to- initial-time coordinate 
which is rewritten as an integration over the initial velocity ttq — dt(po- 

The integral obtained in Eq. (|14ip performs a classical average over a distribution which is 
given by the initial phase space distribution VF[iy9o7 '"'o] propagated to any later time of interest 
by means of the classical equation of motion. For example, assume the path integral to be used 
for the evaluation of the mean value of <P{t) at some time t > to, 



{'P{t)) = / d^od7roM^[¥>o,7ro] J V"p p{t) 5[{d'', + lo'') ^] 

d^odTTo VF[^o, TTo] j P'V ^{t) JW[V - <^ci], (143) 

where tpcX is the solution of the classical equation of motion [df -\-ijj'^)p{t) ~ given initial values 
ifQ and TTo for coordinate and velocity, respectively. The Jacobean J[(p] = |det((5^5[iy9, 5Lp5lp)\ 
is an irrelevant normalisation constant which has been discussed in detail in Ref. [198j and 
references therein. The integral can be performed numerically if W^['/9o: ""o] is a positive definite, 
normalised probability distribution in phase space, in the case discussed here for the initial po- 
sition ipo and momentum ttq of the classical oscillator. Hereby, the functional delta distribution 
implies, that the initial coordinates sampling W are evolved according to the classical oscillator 
equation from to to time t. The resulting distribution is used to perform the average over tp. 

This implies that any correlation function describing the evolution of the quantum me- 
chanical harmonic oscillator can be evaluated using classical statistical sampling given that 
the initial phase-space distribution is positive definite. Nevertheless, the oscillator is of course 
inherently quantum since its ground state does not correspond to the classical ground state 
with Lp = dtp> = 0. A real initial state contains in any case these quantum fluctuation effects, 
though they may be negligible compared to the high-energy classical occupation. Despite this, 
however, no further quantum fluctuations occur during the time evolution of the oscillator. 

The reason for this classicality is, clearly, that the Hamiltonian is quadratic in (p. Consider, 
e.g., the anharmonic oscillator, with classical action 



s[ip] = y" 



(144) 



The quartic term no longer allows to rewrite the action on the close time path into a functional 
linear in (p. On rather obtains 

dt {p+' - =ljdt (^^^3 ^ 1-3^^ _ (145) 

Neglecting the second term which is cubic in ip one can still integrate out ip and obtains the 
classical statistical integral 

d^od7rol^[<^o,7ro] f V"ip 5[{d^ + tu^) ip + Xip% (146) 



again with the classical equation of motion in the argument of the functional delta distribution. 
This procedure is, however, no longer possible if the second term in Eq. p45p is kept, and 
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Fig. 23. (color online) (a) Diagrammatic representation of the correlators in the tp-ip basis. A full 
line indicates the 1- or (p-component, a broken line the 2- or (^-component. Fab{x,y) is the statistical 
correlation function, G^,,{x,y) = Pab{x,y)6{xo - yo) and G'^f,{x,y) = -pab{x,y)0{yo - xq) the retarded 
and advanced Green functions, respectively. Their representation in terms of the real-valued spectral 
correlation function pab{x,y) exposes the 6- functions which imply the respective time ordering in xq, 
yo. (b) Diagrammatic expansion of the quantum vertex contributing to the action as in Eq. (|149p . The 
classical action is lacking the second contribution (red square). 



the explicit path integration over must be kept during the evolution. Hence, the quantum 
anharmonic oscillator indeed shows fluctuation effects beyond classical statistics. 

We finally point out that the classical approximation neglecting the i^^</3 vertex is equivalent 
to the so-called Truncated Wigner Approximation known in quantum optics. See Ref. [199| for 
details as well an approach to include corrections caused by quantum fluctuations analogous to 
what is described above. Closely related results were recently presented in Ref. |200j . 



5.3 Field theory 

We will now summarise the difference between classical and quantum statistical evolution il- 
lustrated in the previous section for the case of fleld theory. Speciflcally, we shall use, in this 
context, the 2PI functional approach introduced in Sect. [3] and study the dynamics of the 
one-dimensional Bose gas. 

Consider the scalar theory of a single-species interacting Bose gas as defined by the classical 
action (|55p . Introducing a linear transformation of the fields as 

this action assumes the form S'[(/5, (p] = Sq[^p^ (p\ + Si■nt['p^ <^], with the free-field-theory action 

i£, (.4 («S) 

and the interaction part 

5'int[<<5, "^1 = -| j (pa{x)(pa{x)tpb{x)tpb{x) ~ ^ J (pa{x)(pa{x)(pb{x)ipb{x), (149) 

in full analogy to the anharmonic oscillator discussed above. 

It is beyond the scope of this article to present the derivation of quantum and classical 
versions of the 2PI effective action and the equations of motion obtained with this, which can 
be found in Ref. [15] ■ In there, also a more careful analysis of the boundary terms neglected above 
is given. Here, we shall present a concise diagrammatic illustration of the differences emerging. 
For this we note that the CTP propagator matrix is transformed, by R from Eq. (|147p . to 

Fab ~T'G^b\ — p ( ^Ib ^Ib \ pT /irn^ 
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where Fab{x,y) is the statistical correlation function. Eq. ([6|), and G^ij{x,y) — pab{x,y)0{xo — 
yo) and G'^ij{x,y) — ~pab{x,y)0{yo — xq) are the retarded and advanced Green functions, 
respectively, which are directly related to the spectral function p. G^ij{x, y), with s,t £ {+, — }, 
denote the full propagator, with xq and yo evaluated on branch s and t of the CTP, respectively. 
Eqs. (|149p and pSOp suggest that the propagators and bare vertices forming the diagrams which 
occur in the effective action when written in the i?-transformed basis can be diagrammatically 
represented as shown in Fig. [221 See Ref. |35] for details of the derivation. 

We consider the possible one-loop as well as higher-loop, closed 2P1 diagrams which can be 
formed from these constituents. Obviously, solid lines only connect to solid lines, and dashed 
lines to dashed ones. Since both vertices have an odd number of solid and dashed legs, no 
diagram up to the two- loop double-bubble diagram, i.e., up to first order in the coupling can 
contain a quantum vertex, depicted as a (red) diamond in Fig. [531 Hence, up to the HFB mean- 
field approximation, the dynamic equations are of classical statistical nature, i.e., apart from 
quantum effects in the initial state, the ensuing dynamics can be simulated by Monte-Carlo 
methods. 



5.3.1 Rescaled fields 



The classical statistical generating functional, in which Sint is lacking the second term in 
Eq. ()149p . exhibits an important reparametrisation property: If the fluctuating fields are rescaled 
according to 

(Pa{x) Lp'^{x) ^ y^(pa{x), >fa{x)'^>p'a{x)^{l/^)'fia{x) (151) 

then the coupling g drops out of S'^^[ip,(p] ~ So[(p,(p] + Sf^^^[(p,0\ where S"^' denotes the term 
in Eq. p49p linear in ip. The free part So[^t0\ remains unchanged and the interaction part 
becomes 

S^L[V',<P'] = -\j^^'AxWa{xWb{xWb{x)- (152) 

Moreover, the functional measure is invariant under the rescaling (|15ip . and the sources can be 
redefined accordingly. Therefore, the classical statistical generating functional becomes inde- 
pendent of g, except for the coupling dependence entering the probability distribution fixing the 
initial conditions. Accordingly, as is well known for the Gross-Pitaevskii equation, the coupling 
does not enter the classical dynamic equations for correlation functions. All the g-dependence 
enters the initial conditions which are required to solve the dynamic equations. 

In contrast to the classical case, this reparametrisation property is absent for the quantum 
theory. After the rescahng ()15ip one is left with S[ip' , (p'] whose coupling dependence is given 
by the interaction part 

5int[<^',^'] = - 2 y ^'a{x)ip',{x)ip',{x)ip',{x) ^ ^- j ^',{xW,{xW,{xW,{x) , (153) 

according to Eq. p49p . Comparing to p52p one observes that the 'quantum' part of the vertex 
encodes all the ^-dependence of the dynamics. 

The comparison of quantum versus classical dynamics becomes particularly transparent 
using the above rescaling. The rescaled macroscopic field and statistical correlation function 
are given by 

'/'Ka;) = V50a(a;), F'^i,{x,y) ^ gFab{x,y), (154) 

while the spectral function pab{x,y) remains unchanged as Eq. (jl50p suggests and is shown in 
Ref. [48] . Similarly, we define for the statistical self-energy entering the dynamic equations (|74l) 



^Ibi.x,y) = gSab{.x,y). 



(155) 
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Fig. 24. (color online) Bubble chains con- 
tributing to the functions I'p and Ip at NLO 
in the 2PI 1/M expansion, see Fig. [23]for def- 
initions. The chains of type (Ql) and (Q2) 
only appear for a quantum system. Type (C) is 
present both in quantum and classical systems. 
These diagrams exhibit that in each term con- 
tributing to the functions I'p and Ip, there is 
at most one loop involving two correlators F' 
or p. In the classical limit, the loops involving 
two p are suppressed compared to those with 
two F'. 



5.3.2 Quantum vs. classical statistical self-energy 



To identify the precise difference between the quantum and the classical time evolution, details 
about the self-energies are required. In the following, we quote the result for the self-energies 
in NLO of the 2PI 1/M expansion introduced in Sect. 13.6731 The self energies ([84|) . expressed 
in the rescaled variables (|154|) and p55p read: 



with A'p{x,y) 



lFix,y) 



with 



A'p{x,y) S^Ap{x,y) 
-^Ap{x,y) A'p{x,y) 



x,y), etc. The functions I'p and 



Kbix.y) 
-\pab{x,y) 

satisfy 



(156) 



VO 



(157) 



Kbix,y)Kb{x,y) - —Pabix,y)pab{x,y) 



np{x,y) = -'F^^{x,y)pab{x,y). 



(158) 



Consider the bubble-chain diagrams occurring within the NLO 1/Af approximation as de- 
picted in Fig. 1101 Connecting the correlators through the respective vertices, in the ip-(p basis 
of Fig. [531 one finds which types of bubble chains appear. The classes of non-vanishing chains 
shown in Fig. 124! confirm the structure of the functions I'p{x,y) and Ip(x,y) which are deter- 
mined by the integral equation ()157|) : In each term contributing to the diagrammatic expansion 
of these functions there is at most one loop containing two F' or two p correlators, the latter 
resulting from either or G^. In addition to this one finds that only the loop containing two 
p correlators goes with the vertex which is present in the quantum case only. Hence, considering 
the classical dynamics, the 
the same reasons, the contribution 
is absent for the classical statistical theory. These arguments can be easily extended to the case 
that (j) and the terms missing in the classical theory be determined 



g^p^ terms are absent together with the 'quantum' vertex. For 
{g^/A)[Ip{x,y) + Pp{x,y)]pab{x,y) to S[t{x,y), Eq. dH, 



5.4 Classicality criterion 



Summarizing, one concludes that all equations ((73)) - ((75)) . (|86)) - ((88)) . and (|156|) - (|158|) remain 
the same in the classical statistical limit except for differing expressions for the statistical 
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components of the self-energy 

^F:,,{x,y)FU{x,y), (159) 

replacing the respective expression in Eqs. ([84|) . The classical statistical self-energies can be 
obtained from the respective quantum ones by dropping two spectral (p-type) components 
compared to two statistical (i^-type) functions. For vanishing (j) where Pf,p ~ this corresponds 
to the result of Ref. [179j . After the rescaling (|15ip the 'quantum' terms can be directly identified 
since they are the only g-dependent terms, which are absent in the classical statistical theory 
according to the above discussion. As a consequence, for the classical dynamics the effects of a 
larger coupling can always be compensated by changing the initial conditions such that Fg, as 
well as (/>Y^, remain constant. This cannot be achieved once quantum corrections are taken into 
account, since they become of increasing importance with growing coupling or reduced initial 
values for F and (j). 

Eq. (|159p describes the differences between quantum and classical statistical equations of 
motion. In turn one can ask under which conditions these differences are negligible. In that case 
the quantum dynamics can be well approximated by classical statistical dynamics. Analysing 
the bubble one finds that a sufficient condition for the suppression of quantum fluctuations 
compared to classical statistical fluctuations is given by [48] 

3 

\Kb{x,y)Fl.^{z,w) \ > -g'^\pab{x,y)pcd{z,w)\. (160) 

This condition is not based on thermal equilibrium assumptions and holds also for far-from- 
equilibrium dynamics. In particular, it is independent of the value of the macroscopic fleld 0. 
It can be applied, of course, also in thermal equilibrium, for which the statistical and spectral 
correlation functions are related by a fluctuation-dissipation relation, see Eq. (|120p in Sect. 14.21 
For large temperatures, k^T ^ uj — fi, one has |F(°'i^'([j, p)|/g ^ \p^°'i\uj,p)\, i.e., the clas- 
sicality condition is fulfilled for all modes whose occupation number ^ F'^'^^'>' {u),p)/ g is much 
larger than 0(1). 

The equivalent statement also holds for nonequilibrium evolutions whenever it is possible 
to define a suitable 'occupation number' from a space-time or energy-momentum dependent 
proportionality between F and p. Away from equilibrium the situation is often considerably 
more complicated. Strictly speaking the condition p60p must be valid at all times and for all 
space points, or momenta in Fourier space, for the classical and the quantum evolution to agree. 
In practice, however, it needs only be fulfilled for time and space averages. In Sect. [53] we will 
demonstrate how quantum evolution can be approximated for not too late times by classical 
statistical dynamics, if the correlation functions satisfy (|160p at initial time. In order to have 
quantum fiuctuations playing a significant role, one needs, according to our above findings, to 
increase the interaction strength g accordingly or change the phase-space structure by changing 
the external trapping potential. For example, in a one-dimensional trap, an effectively strong 
coupling and strong quantum fiuctuations can be induced by reducing the line density of atoms 
while their interaction strength is kept constant. 



5.5 Quantum vs. classical evolution of an ultracold Bose gas 

In this section we apply the theoretical methods summarised above to the equilibrating ID Bose 
gas studied already in Sects. 14. 3[ and l4.5l Along the lines discussed in the previous sections 
we compare the evolution involving only classical statistical fluctuations with that which also 
takes into account quantum corrections. Both, the classical and the quantum gas are assumed 
to be initially characterised by the same far-from-equilibrium initial conditions far. We solve 
Eqs. ([74]). with the self-energies now given by Eqs. (|159p . for the initial conditions given in 



^abi^^y) 

Yjf / \ classical limit 

llp[x^y) > 



Ultracold gases far from equilibrium 



53 




0.0001 



0.001 0.01 
f[s] 



Fig. 25. (Color online) The nor- 
malised momentum-mode occupa- 
tion numbers n{t; p) /n\L for the clas- 
sical gas (black solid lines) compared 
to their quantum counterparts from 
Fig. [15] (red dashed lines), as func- 
tions of time. Shown are the popu- 
lations of the modes with p = pi = 
2Ns/Lsm{iTT/Ns), i = 0,1, Ns/2, 
and one has n{t;—p) = n(t;p). In 
contrast to the quantum statistical 
evolution there is no quick dephasing 
in the classical case, such that the ini- 
tially empty modes become only sub- 
sequently filled, as is discussed in the 
main text. At large times, the clas- 
sical and quantum gases necessarily 
evolve to different distributions. 



Section [3.8.11 first with a = 0. For comparison we also consider the case a = 1 such that the 
classical and quantum initial correlation functions F and p are identical. 

The time evolution of the initially Gaussian far-froni-cquilibrium momentum distribution 
of the classical gas is shown in Fig. [551 The mode occupations are shown as (black) solid 
lines, and for better comparison, we have quoted the quantum evolution from Fig. [TSjas (red) 
dashed curves. One observes that the time evolution of the modes with occupation number 
n(t;p) > 1, i.e., n{t;p)/niL > 10""^, is, for most of the time, identical to that obtained in 
the quantum case, confirming condition (|160p . As expected, for the strongly populated modes 
of a weakly interacting bosonic gas quantum fluctuations do not play a significant role for 
not too large times. Only when the evolution approaches the equilibrium state, the differences 
between quantum and classical statistics are expected to lead to a Bose-Einstein and classical 
distribution, respectively. We point out that, although we have chosen, for our comparison, 
the same initial occupation numbers in the two cases, the total energies are different since the 
initial correlation functions differ according to Eq. (|108p . Hence, also the final-state occupation 
numbers of the low momentum modes can differ. 

We finally point to the substantial differences in the short-time evolution of the weakly 
populated modes. While, in the quantum gas, the large-momentum mode populations are all 
growing at the same rate, the modes of the classical gas become populated much more gradually. 
The quantum-gas modes are occupied by 0.01 and 1 particle per mode already between 0.5 and 
1 fis, while the classical modes need up to ten times longer. The distinct quantum behaviour of 
the modes can be understood as follows. At the energies present in such an ultracold gas, the 
atomic interactions are essentially pointlike, i.e., the range of the potential is not resolved and 
the coupling function or scattering amplitude is constant over the range of relevant momenta. 
Hence, in a single scattering event, the distance of two atoms is localised to zero, such that the 
relative momentum of the atoms is completely unknown immediately after the collision. This 
means that the transfer probability of the atoms is the same for any final momentum mode, and 
this is observed as a quick, collective population in the respective, so far essentially unoccupied 
modes. 

For comparison we repeated our calculations for a non-local interaction potential. In the 
momentum domain, this corresponds to a coupling function which is cut off at large momenta, 
and we chose the cutoff within the range of the momenta shown explicitly in Fig. 1251 In this 
case we find that the quantum evolution is modified such that it becomes similar to that of the 
classical gas. In particular, all modes above the cutoff populate gradually one after each other. 

The differences between the quantum and classical evolutions shown above depend, however, 
considerably on the choice of initial conditions. To compare the characteristics of the evolutions 
which are independent of the initial choice of F, we have repeated the classical calculations for 
an initial momentum distribution where, as compared to before, a constant occupation number 
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Fig. 26. (Color online) The nor- 
malised momentum-mode occupa- 
tion numbers n{t; p) /niL for the clas- 
sical gas (black solid lines) compared 
to their quantum counterparts from 
Fig. [15] (red dotted lines), as func- 
tions of time. All parameters are cho- 
sen as in Fig. 1251 except for a = 
6.5 ■ 10" m"\ The dashed (blue) lines 
show n{t,p)~l/2, i.e., after subtract- 
ing the additional flat initial distri- 
bution which simulates the quantum 
"zero point fluctuations". 

1/2 has been added. Hence, we chose a = 1 in the initial values of F, Eq. (|108|) . as in the 
quantum case, such that F is identical for xq — yo ^ io in the classical and quantum cases. The 
results are shown in Fig.[2Sl The (red) dotted lines show, again, the quantum evolution, while the 
classical mode populations for the same initial conditions for F and p are shown as solid (black) 
lines. Subtracting 1/2 from each Faa{t,t;p) gives the dashed (blue) lines. One finds, that during 
the initial period the evolution of the variation of the high-momentum modes with respect to 
their initial occupation is identical to the quantum evolution of the occupation numbers. At 
intermediate times, however, there are deviations which lead to occupation numbers up to 
a percent lower than 1/2. Although the chosen initial conditions which, in the classical case 
correspond to a base occupation of each mode with 1/2 atom, seem unphysical, our results 
show that for the dilute, weakly interacting gas under consideration, there are differences only 
in those modes which, in the mean, are populated with less than one atom. Quantum statistical 
fluctuations play a role only for these modes. 



6 Summary 

We have discussed functional field theoretical approaches to far-from-equilibrium quantum 
many-body dynamics. The central topic were methods based on the real-time two-particle 
irreducible (2PI) effective action to derive coupled integro-differential equations of motion for 
correlation functions. These equations were used to study, as a working example, the long- 
time evolution and equilibration of an ultracold Bose gas in one spatial dimension, starting 
in a state far away from equilibrium. They allow to access phenomena which are beyond the 
reach of transport or kinetic theory, i.e., beyond that of (quantum) Boltzmann equations. Non- 
perturbative approximation schemes were discussed which allow to gain access to the dynamics 
of strongly interacting systems and to long-time evolution and thermalisation. 

By use of functional renormalisation-group (RG) techniques the 2PI dynamic equations 
were rederived in a particular truncation scheme of the RG flow equations, and an alternative 
non-perturbative approximation for use with strongly correlated systems was obtained. 

The transition from initial to late-time nonequilibrium evolution and the emergence of trans- 
port theory were studied in detail and illustrated with the one-dimensional Bose gas example. 
In the last part, we worked out the distinction between classical and quantum statistical fluc- 
tuations during the time evolution of a quantum system, first on the example of a quantum 
mechanical harmonic oscillator described by a Feynman path integral. Extending the formula- 
tion to functional integrals allowed us to identify the effect of quantum fluctuations in terms of 
vertex contributions which are absent in the classical case. 

Nonequilibrium quantum field theory has seen a vigorous development during the past 
decade, but many questions remain open leaving space for further progress. Nonperturbative 
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approximations to describe strongly correlated systems are at hand but need to be explored in 
more detail concerning their applicabihty to various specific configurations. For such applica- 
tions also numerical methods are desirable, as they are existing and used for classical statistical 
evolutions. New field theoretical methods which provide the nonequilibrium time evolution, 
e.g., of equal-time correlation functions but do not require memory integrals over the history 
of the non-equal-time functions to be evaluated explicitly would help towards applications. 

Experiments with ultracold gases have the potential to provide precise information about 
nonequilibrium dynamics of correlations. The methods described in this article must be made 
accessible to these experiments, i.e., concrete applications need to be developed and optimal 
observables need to be found. With these tasks in mind dynamical field theory is expected to 
continue to be and develop even stronger to an exciting research field. 
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